arXiv: 1506.0103lv2 [astro-ph.CO] 13 Oct 2015 


Mon. Not. R. Astron. Soc. 000,[Tp4l(2015) Printed 15 October 2015 (MN MeX style file v2.2) 


Towards the statistical detection of the warm-hot intergalactic 
medium in inter-cluster filaments of the cosmic web* 

Nicolas Tejos/f J. Xavier Prochaska/ Neil H. M. Crighton,^ Simon L. Morris,^ 
Jessica K. Werk/ Tom Theuns,^’^ Nelson Padilla,^ Rich M. Bielby^ and Charles W. 
Finn^’^ 

^ Department of Astronomy and Astrophysics, UCO/Lick Observatory, University of California, 1156 High Street, Santa Cruz, CA 95064, USA 
^ Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Hawthorn, Victoria 3122, Australia 
® Department of Physics, Durham University, South Road, Durham, DHl 3LE, UK 

^ Institute for Computational Cosmology, Department of Physics, University of Durham, South Road, Durham, DHl 3LE, UK 
® Instituto de Astrofi'sica, Centro de Astro-lngenien'a, Pontificia Universidad Catolica de Chile, Av. Vicuna Mackenna 4860, Santiago, Chile 


Accepted version, 15 October 2015 


ABSTRACT 

Modem analyses of structure formation predict a universe tangled in a ‘cosmic 
web’ of dark matter and diffuse baryons. These theories further predict that at low-z, a 
significant fraction of the baryons will be shock-heated to T ^ 10® — 10^ K yielding 
a warm-hot intergalactic medium (WHIM), but whose actual existence has eluded a 
hrm observational conhrmation. We present a novel experiment to detect the WHIM, by 
targeting the putative hlaments connecting galaxy clusters. We use HST/COS to observe 
a remarkable QSO sightline that passes within Ad = 3 Mpc from the 7 inter-cluster axes 
connecting 7 independent cluster-pairs at redshifts 0.1 ^ z ^ 0.5. We find tentative 
excesses of total H I, narrow H I (NLA; Doppler parameters b < 50 km s“^), broad 
H I (BLA; b ^ 50 km s and O VI absorption lines within rest-frame velocities of 
Av < 1000km s“^ from the cluster-pairs redshifts, corresponding to ^ 2, ^ 1.7, ^ 6 
and ^ 4 times their field expectations, respectively. Although the excess of O VI likely 
comes from gas close to individual galaxies, we conclude that most of the excesses of 
NLAs and BLAs are truly intergalactic. We find the covering fractions, fc, of BLAs 
close to cluster-pairs are ~ 4 — 7 times higher than the random expectation (at the ~ 2cr 
C.I.), whereas the of NLAs and O VI are not significantly enhanced. We argue that a 
larger relative excess of BLAs compared to those of NLAs close to cluster-pairs may be 
a signature of the WHIM in inter-cluster hlaments. By extending the present analysis to 
tens of sightlines our experiment offers a promising route to detect the WHIM. 

Key words: -intergalactic medium -quasars: absorption lines -large scale structure of 
the Universe -galaxies: formation 


1 INTRODUCTION 

Perhaps the most distinctive feature of the cosmic web is its 
intricate pattern of filamentary structures. Cosmological simu¬ 
lations in a A cold dark matter (ACDM) paradigm predict that 
these filaments account for ~ 40% of all mass in the Universe at 

2 = 0 and occupy roughly ~ 10% of the volume (e.g.|Arag6n-| 
|Calvo et al.|2010t . When gas and hydrodynamical effects are 
included in these simulations, a remarkable conclusion is reached: 
~ 30 — 50% of baryons at low -2 should reside in dense filaments, 
primarily in the form of a diffuse gas phase with temperatures 
T ~ 10® — 10^ K, which would be very difficult to detect (e.g. 
|Cen & Ostriker|19^^|Dave et al.|20OT| l. This material is usu¬ 
ally referred to as the warm-hot intergalactic medium (WHIM), 


* Based partly on observations made with the NASA/ESA Hubble Space 
Telescope under program GO 12958. 
f E-mail: ntejos@gmail.com 


and indeed is currently the best candidate to host a significant 


fraction of the so-called ‘missing baryons’ at 2 < 1 (Persic & 


Salucci|19^JFukugitae£ayi9^ Bregman|2007[[Prochaska 

& Tumlinson|2009| [Shull et al.|2012| and references therein). 

According to these models, the physical origin of the WHIM is 
through gravitational shocks from the collapse of matter into the 
large scale structure (LSS) of the Universe. 

One well-studied example of gravitational shock-heating 
is the so-called intracluster medium (ICM) of galaxy clusters, 
where the virial temperatures typically reach T ~ 10^ — 
10* K. A plasma at these temperatures mostly cools through 
Bremsstrahlung (a.k.a. free-free) thermal radiation, emitting X- 
rays at ~ keV energies that may be observed with modem satel¬ 
lites (e.g. IKravtsov & Borgani||2012] and references therein). 
X-ray spectroscopy has also revealed the presence of high- 
ionization state metal emission lines in the ICM, consistent with 
these large temperatures (e.g. [Sanders et al.|2008^ . Thereby one 
constrains the density, chemical abundances and morphology 
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of the ICM. Several decades of research have revealed a highly 
enriched medium (~ | solar) with a total mass consistent with 
the cosmic ratio of baryons to dark matter (e.g. |Allen et al.|2008^ . 

In the ACDM paradigm, galaxy clusters correspond to the 
nodes of the cosmic web, i.e., they mark the intersection of sev¬ 
eral filamentary threads. These models further predict that matter 
flows through the filamentary structures, driving the growth of 
the galaxy clusters. Ideally, one would image these filaments in a 
similar manner to the ICM to reveal their structure and physical 
properties as tests of the cosmic web paradigm. Unfortunately, 
once at the outskirts of galaxy clusters, the densities and temper¬ 
atures are too low for viable X-ray detection in emission (e.g. 
Bremsstrahlung radiation is proportional to the density squared 
of the emitting gas). To study this dominant component of the 
cosmic web and its putative relationship to a WHIM, one must 
pursue alternate strategies. 

In principle, one may scour the volumes surrounding galaxy 
clusters for signatures of cosmic filaments. A random search, 
however, would be compromised by the fact that their volume 
filling factor is predicted to be low, even in this environment. To 
raise the probability of isolating a cosmic filament, researchers 
have turned to pairs of neighbouring clusters on the expectation 
that these massive structures will be preferentially connected. 
Indeed, cosmological dark matter simulations find high proba¬ 
bilities of having a coherent filamentary structure between close 
(< 20Mpc) and massive (> 10^"^ Mq) galaxy clusters (e.g. 


Colberg et al.|2005| [Gonzalez & Padilla|20T^ |Arag6n-Calvo] 

et al. 12010 1 . This probability is mostly a function of the galaxy 
cluster masses and the separation between them: the larger the 
masses and the shorter the separation, the higher the probability. 
Therefore, the volume between close pairs of galaxy clusters 
is a natural place to search for signatures of filaments and an 
associated WHIM. 

Inter-cluster filaments (i.e. filaments between galaxy cluster 
pairs) have been inferred from galaxy distributions, either indi¬ 
vidually from spectroscopic galaxy surveys (e.g. |Pimbblet et ak] 
|2004^, or by stacking a nalysis from photometric galaxy surveys 
(e.g. Zhang et al.|2013^ . While these studies confirm the strategy 
to focus on cluster pairs, they provide limited information into 
the nature of cosmic filaments; these luminous systems represent 
< 10% of the baryonic matter, their distribution and motions 
need not trace the majority of the gas, and they offer no insight 
into the presence of a WHIM. 

Promising results from stacking multiple inter-cluster 
regions have found an excess of X-ray counts in such regions 
with respect to the background jFraser-McKelvie et al.|20TT) . 
In contrast to galaxies, one would be truly observing the 
bulk of baryonic matter. Unfortunately, the geometry of the 
emission and the actual origin of the detected photons was not 
well constrained by this original work. Remarkable detection 
of individual inter-cluster filaments have also been reported 


from gravitational weak lensing signal i Dietrich et al. 


and X-ray emission ( |Kull & Bdhring^ 1999[ Werner et al.| 
|2008| l. Despite their indisputable potential for characterizing 
cosmological filaments, these techniques are currently limited 
to the most massive systems with geometries maximizing the 
observed surface densities, i.e. filaments almost aligned with the 
line-of-sight (LOS). 

To complement these and other relevant studies to address 
the ‘missing baryons’ problem (e.g. [Nevalainen et al.||2015[ 
[Planck Collaboration et al.j2015[pHemandez-Monteagudo et ^ 


|2015| l, we have designed a program to detect the putative fila¬ 
ments connecting cluster pairs in absorption. This technique has 
several advantages over attempts to detect the gas in emission. 
First, absorption-line spectroscopy is linearly proportional to the 
density of the absorbing gas, offering much greater sensitivity 
to a diffuse medium. Second, the absorption lines encode the 
kinematic characteristics of the gas, including constraints on 
the temperature, turbulence, and LOS velocity. Third, one may 
assess the chemical enrichment and ionization state of the gas 
through the analysis of multiple ions. The obvious drawback to 
this technique is that one requires the fortuitous alignment of a 
bright background source with these rare cluster pairs, to probe 
a greatly reduced spatial volume: in essence a single pinprick 
through a given filament. However, with a large enough survey 
one may also statistically map the geometry/morphology of the 
filaments. 

Here we focus on far ultra-violet (FUV) spectroscopy lever¬ 
aging the unprecedented sensitivity of the Cosmic Origins Spec¬ 
trograph (COS) onboard the Hubble Space Telescope (HST), to 
greatly increase the sample of inter-cluster filaments probed]^ 
With such UV capabilities we can directly access H I Lya— 
the strongest and most common transition for probing the 
intergalactic medium (IGM). Having direct coverage of H I inde¬ 
pendent of the presence of metals is of great value for detecting 
the WHIM (e.g. [Richter et al.|2006[[Danforth et al.|20I0^ , be¬ 
cause this medium may remain metal poor. Neutral hydrogen 
generally traces cool and photoionized gas, but it may also trace 
collisionally ionized gas in the WHIM through broad (Doppler 


parameters b > 50 km s ^) lines (e.g. Tepper-Garcfa et al. 


2012 '. 


Although the circumgalactic medium (CGM) surrounding galax¬ 
ies is responsible for producing HI absorption lines (especially at 


column densities > 10^® cm 


e.g. 


Tumlinson et 


akp^, 

majority of them must arise in the diffuse IGM (e.g. jProc haskaj 
jet aI.|20lT]|Tejos et al.|20I2||20I4^ . FUV spectroscopy also al¬ 
lows the detection of the O VI doublet, a common highly ionized 
species. The physical origin of O VI absorption lines is contro¬ 
versial, including scenarios of photoionized and/or collisionally 
ionized gas in the CGM of individual galaxies and/or galaxy 
groups (e.g. Tripp et al.|2008[ Thom & Chen|2008 Wakker 


the 


|Savage|2009[ Stocke et al.|2014[ S^age et al.|2014) . Thus, a 

collisionally ionized component could well be present, some of 
which may come from a WHIM (although seejOppenheimer (&| 
|Dave|20^ [Tepper-Garcfa et al.|2Ql l) . 

In a more general context, H I and O VI offer an optimal 
approach to study filamentary gas in absorption. As mentioned, 
this pair of diagnostics correspond to the most common tran¬ 
sitions observed in the low-z Universe (e.g. Danforth & Shullj 


[2008[rinpp et al.|2008t[Danforth et al.|2014 i, allowing a good 

characterization of the background signal against which one may 
search for signatures of WHIM in filamentary gas. Such signa¬ 
tures could include an elevated/suppressed incidence, covering 
fractions, and/or unique distributions in the strengths or widths 
of the absorption features. In contrast to studies where absorption 
syste ms could be associated with filaments on an individual basis 
(e.g. Aracil et al.|2006 Narayanan et al. 2010|lp]our method- 


^ We note that X-ray spectroscopy could also be used to trace the 
WHIM in absorption, mostly through 0 VII absoiption lineMe^gJNicasttol 


et al. 2005[[Fang et al.|2010[ [Zappacosta et al.|2010[ [Nicastro et ^ 


^ See also [jauzac et al.[p012j for a weak lensing signal of a filament 
connecting to a single galaxy cluster. 


2010 . However, the poor sensitivities of current X-ray spectrographs 
considerably limits the sample sizes for these studies. Furthermore, such 
poor sensitivities and poor spectral resolutions make the interpretation of 
signals particularly challenging (e.g.[Yao et al.|2012^. 

® We note that even in individual cases where such absorption does 
coincide with known structures traced by galaxies, it is still unclear 
whether the gas is actually produced by a WHIM or individual galaxy 
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Table 1. Properties of the observed quasi-stellar object (QSO) Q1410. 


QSO Name 

(1) 

R.A. 

(hr min sec) 
(2) 

Dec. 

(deg min sec) 

(3) 

2:qso 

(4) 

r 

(5) 

Magnitudes 

NUV 

(6) 

FUV 

(7) 

SDSS J141038.39-I-230447.1 

14 10 38.39 

+23 04 47.18 

0.7958 

17.C 

1 17.4 

18.7 


(1) Name of the QSO. (2) Right ascension (J2000). (3) Declination (J2000). (4) Redshift of the 
QSO. (5) Apparent r (visual) magnitude from Sloan Digital Sky Survey (SDSS). (6) Apparent 
near-ultra-violet (UV) magnitude from GALEX. (7) Apparent far-UV magnitude from GALEX. 


ology is statistical in nature and a large sample of independent 
structures must be collected. 

The current advent of big extragalactic surveys makes our 
approach feasible. For instance, the SDSS ( |Ahn et ST]|2014| | 
provides large samples of LSS traced by galaxies and known 
QSOs in the same volume. In particular, by using the galaxy 
cluster catalog of |Rykoff et al.|f2014^ we have constructed a 
cluster-pair sample and found that, on average, a random sightline 
extending between 0.1 Si 2 Si 0.5 intersects 1 ± 1 independent 
cluster-pairs with projected separations of Si 3 Mpc to the inter¬ 
cluster axis (defined as the line segment joining the centers of 
the two galaxy clusters of a pair; see Figure[T]for an illustration, 
and Sectionj^for further details), with a very skewed distribution 
towards zero (see Section[3^. In order to enhance the efficiency 
however, we have cross-matched such cluster-pair sample with 
known FUV luminous QSOs from the |Schneider et aT] ( |2010[ l 
catalog, and identified particular sightlines intersecting more than 
one of these structures. Our approach is highly complementary 
to that of |Wakker et al.|j2015^ , where a single galaxy filament is 
targeted with multiple QSO sightlines. 

In this paper we present HST/COS FUV observa¬ 
tions of a single bright QSO at 2 « 0.8 (namely SDSS 
J141038.39-1-230447.1, hereafter referred to as Q1410; see 
Table [TJ, whose unique sightline intersects 7 independent 
cluster-pairs within 3 Mpc from their inter-cluster axes. This 
sightline is highly exceptional; the random expectation of finding 
such a number of cluster-pairs is < 0.01% (see Section 
With this one dataset, we offer a first statistical assessment of the 
presence of diffuse gas close to cluster-pairs. Although we are 
only reporting tentative results (~ 1 — 2cr c.l.), the primary focus 
of this manuscript is to establish the experimental design and 
methodology. Future work will extend the present study to tens 
of sightlines, eventually leading towards the statistical detection 
of the WHIM in inter-cluster filaments of the cosmic web. 

Our paper is structured as follows. In Section|^we present 
both the galaxy cluster catalog used to create our cluster-pair 
sample and our HST/COS observations of Q1410. In Sectionj^ 
we characterize the volume around the Q1410 sightline in terms 
of known clusters and cluster-pairs, quantifying how unusual the 
Q1410 field is. In Section|^we provide a full characterization 
of the Q1410 HST/COS FUV spectrum in terms of absorption 
line systems, regardless of the presence of known intervening 
structures. In Sectionj^we present our methodology to cross¬ 
match the information provided by the cluster-pair sample and 
absorption line systems, while in Section|^we present our ob¬ 
servational results for the Q1410 field. A discussion of these 
results is presented in Section]^ and a summary of the paper 
is presented in Section]^ Supplementary material is presented 


halos (e.g.|Stocke et al.|2006||Prochaska et al.|2011||Tumlinson et al.| 
[M7|[WlIiams et al.|2013y~~^ 


in the Appendix. All distances are in co-moving coordinates as¬ 
suming Ho = 67.3km s“^ Mpc“^, Dm = 0.315, Qa = 0.685, 
fe = 0 (unless otherwise stated), where Hq, Dm, Da and k are 
the Hubble constant, mass energy density, ‘dark energy’ density 
and spatial curvature, respectively ( [Planck C ollaboration et al.| 
[20T4t . 


2 DATA 

2.1 Galaxy clusters 

In this section we briefly describe the cluster catalog used in 
the present paper. We used red-sequence Matched-filter Proba¬ 
bilistic Percolation (redMaPPer) ( jRykoff et al.|2014 |l appl ied to 
the SDSS Eighth Data Release (DR8) ( Aihara et al.|201l| (. This 
is one of the largest galaxy cluster catalogs currently available, 
containing ~ 25 000 ricl^alaxy clusters (> 20 galaxies having 
luminosities L ^ 0.2L*|jat 0.08 ^ 2 Si 0.55. 

The redMaPPer catalog is very well suited for statistical 
analysis: it defines clusters properties in terms of probabilities 
(e.g. position, richness, redshift, galaxy members), with a well 
understood selection function; it adopts an optimal mass-richness 
relationship ( jRykoff et al.|201^ ; and it has high completeness 
and purity levels compared to others cluster catalogs jRykoff] 
jet al.|2014[rRozo & Rykoff|2014[|Rozo et al.|2014| l. 

In this paper we used an extension of the published 
redMaPPer catalog, including galaxy clusters with richness be¬ 
low 20 but larger than loj^The mass-richness relation relevant 
to the redMaPPer catalog is. 


In 


M200 


= 1.72 + 1.08 In ( — 


( 1 ) 


hjQ lO^^Af0 

with a typical scatter of ~ 0.25 in ln(M) ( Rykoff et al.|2012 1 , 
where M 200 is the total mass enclosed within an overdensity of 
200 times the critical density of the Universe,pcrit; /170 is the di¬ 
mensionless Hubble parameter hro = Ho/ [70 km s“^ Mpc“^); 
Mq is the mass of the sun; and A is the richness of galaxies 
with luminosities L 0.2L* (corrected for incompleteness). 
Extrapolating this relation to A = 10 we get a minimum mass 
A /200 (A = 10) ~ 0.8 X 10^'^Mq in our assumed cosmology 
(see end of Section[T]l. Therefore, our adopted limit should still 
ensure a reasonable large minimum mass limit. We also com¬ 
puted an estimate of the virial radius, R 200 , simply defined as 
the radius at which the M 200 is enclosed. 


R 200 = 


M21 


7r200pc: 


1/3 


( 2 ) 


^ According to their calibration, this richness limit correspo nds to a mass 
of M ~ 1.8 10^^ Mq (uncertain up to ~ 0.25 in ln(A7); Rykoff et al. 
|2012| see also Section [+T| . 

^ Kindly provided by E. Rykoff and E. Rozo (private communication). 
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Projected in the sky: 


QSO 


Cluster 




Inter-cluster axis 


Cluster 

-o 


Along the line-of-sight: 


Cluster 

o 


Inter-cluster axis 



Figure 1. Schematic representation of our survey geometry projected in the sky (top) and along the line-of-sight (bottom), for a single cluster-pair. 
Galaxy clusters are represented by red circles, while the inter-cluster axis is represented by the grey dashed lines. The QSO itself is represented by a 
eight-pointed star. The impact parameter between the QSO sightline and the inter-cluster axis is defined as Ad, while the distance to the closest cluster of 
the pair along the projected inter-cluster axis is defined as Aa;. A rest-frame velocity window around the position of the cluster-pair is defined as it Au. 


We note that because we are only using these clusters as 
tracers of high density regions in the cosmic web, the exact mass 
of the clusters will not be particularly relevant, making potential 
systematic uncertainties in the mass-richness calibration and its 
extrapolation to lower values not a critical issue. Moreover, low 
richness cluster samples suffer more from impurity and incom¬ 
pleteness ( [Rykoff et al.|2014k but such issues should not create a 
fake signal in differentiating the properties of absorption line sys¬ 
tems in environments traced by these cluster-pairs with respect 
to the field, in the presence of real inter-cluster filaments. On the 
contrary, our approach is conservative in the sense that impurity 
and incompleteness would dilute such a signal, if present. 

For the purposes of this paper, we also required the clusters 
to lie between 0.1 ^ ^ 0.5, yielding a total of 162 144 clusters 

(see Section [34l and Appendix|^for further details). 

Although the redMaPPer catalog is mostly based on pho¬ 
tometry, the spectroscopic redshift of the most likely centre is 
also given when available (typically from its brightest cluster 
galaxy (BCG))|^This is advantageous for our experiment; a high 
precision in the cluster redshifts is needed for a reliable associa¬ 
tion with the gas observed in absorption with inter-cluster fila¬ 
ments (if any). About 44% (71 701/162 144) of the redMaPPer 
clusters have spectroscopic redshifts, yielding velocity precision 
of ~ 30 km s“^ in those clusters’ rest frames. We note that for 
the subsample of clusters most relevant to our paper (i.e. within 


® We note that the typical photometric redshift uncertainties for the 
redMaPPer clusters are of the order of 5z ~ 0.006. 


20Mpc of the Q1410), such fraction increases to ~ 70% (40/57; 
see Section [jT] for further details). 

2.2 Q1410 absorption lines 

2.2.1 Specific selection of Q1410 

In this section we describe in detail the selection criteria of our 
targeted QSO, Q1410. We emphasize that the original selection 
of Q1410 was done using the galaxy cluster catalog published by 
Hao et al.| ( |2010^ instead of the redMaPPer catalog jRykoff et al.| 
20 14^ used here. 

The Gaussian Mixture Brightest Cluster Galaxy (GMBCG) 
catalog ( |Hao et al.|2010^ is based on data from the Seventh Data 
Release (DR7) of the SDSS ( [Abazajian et al.|2009^ , and simi¬ 
lar to redMaPPer, it is mainly based on photometry. From the 
GMBCG catalog we searched for cluster-pairs where at least 
one member has a spectroscopic redshift and where the redshift 
difference between them is less than 3x the combined redshift 
uncertainty, and over clusters having GMBCG richness > 15. 
We then measured the transverse co-moving separation between 
clusters at the redshift of the cluster with spectroscopic identifi¬ 
cation (if both clusters had spectroscopic redshifts we used the 
average redshift), and kept the ones separated by < 25 Mpc. 

We selected our target from the QSO catalog published 
by |Schneider et aL] ( |2010[ l, which is also based on SDSS DR7 
data. This catalog comprises > 100 000 QSOs with well known 
magnitudes and spectroscopic redshifts. We looked for QSOs 
having redshifts greater than individual cluster-pairs and located 
inside their sky-projected cylinder areas as defined above. We 


© 2015 RAS, MNRAS 0Q0,[T|[24l 





















Towards detecting the WHIM in inter-cluster filaments 


5 


Table 2. Properties of the redMaPPer clusters at 0.1 ^ ^ sg 0.5 and within 20Mpc from the Q1410 sightline. 


Cluster ID 

CD 

R.A. 

(degrees) 

(2) 

Dec. 

(degrees) 

(3) 

2^spec 

(4) 

•^photo 

(5) 

Richness 

(6) 

Mass 

(ioi'‘ Mq) 
(7) 

-^200 

(Mpc) 

(8) 

Impact parameter 
(degrees) (Mpc) 

(9) (10) 

C-R200) 

(ID 

1 

212.994666 

21.418861 

0.1335 

0.150 ± 0.006 

11.0 

0.92 

0.9 

1.6895 

16.9 

18.5 

2 

213.481691 

22.808351 

0.1376 

0.147 ± 0.005 

20.4 

1.80 

l.I 

0.8039 

8.3 

7.3 

3 

214.114059 

23.256258 

0.1381 

0.138 ± 0.005 

33.0 

3.03 

1.4 

1.3484 

13.9 

10.3 

4 

213.817562 

24.020988 

0.1386 

0.152 ± 0.006 

16.5 

1.43 

1.1 

1.4184 

14.7 

13.9 

5 

213.014763 

22.125043 

0.1413 

0.141 ± 0.006 

18.3 

1.60 

1.1 

1.0093 

10.6 

9.7 

6 

213.223623 

22.319874 

0.1417 

0.150 ± 0.005 

20.3 

1.79 

1.1 

0.9208 

9.7 

8.6 

7 

213.330826 

22.510559 


0.143 ± 0.005 

17.9 

1.56 

1.1 

0.8405 

9.0 

8.3 

8 

213.664909 

22.216942 

0.1535 

0.171 ± 0.007 

10.9 

0.91 

0.9 

1.2667 

14.5 

16.0 

9 

212.712357 

23.086999 

0.1580 

0.157 ± 0.006 

13.0 

1.11 

1.0 

0.0487 

0.6 

0.6 

10 

212.162442 

21.946327 

0.1596 

0.145 ± 0.006 

11.3 

0.95 

0.9 

1.2231 

14.5 

15.9 

11 

212.999309 

24.435343 

0.1612 

0.171 ± 0.006 

16.1 

1.39 

1.0 

1.3907 

16.6 

16.1 

12 

212.674015 

22.495683 

0.1733 

0.179 ± 0.006 

23.0 

2.05 

1.2 

0.5842 

7.5 

6.4 

13 

212.895979 

23.605206 

0.1787 

0.167 ± 0.007 

11.4 

0.96 

0.9 

0.5684 

7.5 

8.3 

14 

213.557346 

22.914918 

0.1918 

0.195 ± 0.007 

12.4 

1.05 

0.9 

0.8423 

11.9 

12.8 

15 

212.123899 

22.560692 

0.2220 

0.241 ± 0.012 

11.2 

0.94 

0.9 

0.7167 

11.6 

13.1 

16 

211.805856 

23.845513 

0.2359 

0.233 ± 0.008 

21.4 

1.90 

l.I 

1.0955 

18.8 

16.8 

17 

213.277561 

24.079481 

0.2392 

0.234 ± 0.012 

10.6 

0.89 

0.9 

1.1488 

20.0 

23.1 

18 

213.029395 

23.923686 

0.2424 

0.254 ± 0.010 

24.9 

2.23 

1.2 

0.9094 

16.0 

13.6 

19 

212.319446 

22.314272 

0.2914 

0.482 ± 0.021 

11.4 

0.96 

0.9 

0.8275 

17.3 

19.8 

20 

212.415280 

22.627699 

0.3127 

0.306 ± 0.020 

11.4 

0.96 

0.9 

0.5052 

11.3 

13.0 

21 

213.395595 

23.114057 


0.327 ± 0.020 

12.2 

1.03 

0.9 

0.6775 

15.8 

17.9 

22 

212.831535 

22.621344 

0.3382 

0.368 ± 0.017 

30.0 

2.73 

1.2 

0.4849 

11.6 

9.6 

23 

212.723980 

22.929555 

0.3412 

0.357 ± 0.017 

28.1 

2.54 

1.2 

0.1614 

3.9 

3.3 

24 

213.000928 

22.437411 


0.345 ± 0.018 

23.6 

2.10 

l.I 

0.7152 

17.5 

15.7 

25 

212.293470 

23.265859 

0.3465 

0.349 ± 0.019 

21.2 

1.88 

l.I 

0.3849 

9.4 

8.8 

26 

212.684410 

23.402739 

0.3506 

0.358 ± 0.024 

11.4 

0.96 

0.9 

0.3237 

8.0 

9.4 

27 

212.770942 

23.232412 

0.3508 

0.357 ± 0.022 

10.4 

0.87 

0.8 

0.1836 

4.6 

5.5 

28 

213.066601 

23.075991 

0.3511 

0.378 ± 0.021 

20.6 

1.82 

1.1 

0.3741 

9.3 

8.8 

29 

212.600828 

23.822434 

0.3520 

0.358 ± 0.020 

15.0 

1.29 

0.9 

0.7446 

18.5 

19.7 

30 

212.854438 

22.531418 


0.365 ± 0.021 

13.5 

1.15 

0.9 

0.5769 

14.8 

16.5 

31 

212.861147 

22.341353 

0.3718 

0.385 ± 0.016 

22.8 

2.03 

1.1 

0.7614 

19.9 

18.3 

32 

213.011303 

23.110255 

0.3722 

0.395 ± 0.019 

13.6 

1.16 

0.9 

0.3246 

8.5 

9.4 

33 

212.737954 

23.677308 

0.3725 

0.379 ± 0.017 

18.4 

1.61 

1.0 

0.6018 

15.8 

15.7 

34 

213.115594 

22.812719 

0.3725 

0.400 ± 0.022 

16.9 

1.46 

1.0 

0.4973 

13.0 

13.4 

35 

212.604626 

23.277634 


0.376 ± 0.018 

22.2 

1.98 

1.1 

0.2043 

5.4 

5.0 

36 

212.538462 

23.376505 


0.384 ± 0.019 

17.6 

1.53 

1.0 

0.3170 

8.5 

8.7 

37 

212.800638 

23.052477 

0.4138 

0.426 ± 0.015 

20.0 

1.76 

1.0 

0.1323 

3.8 

3.7 

38 

212.609281 

22.997888 

0.4159 

0.425 ± 0.015 

10.6 

0.89 

0.8 

0.0942 

2.7 

3.4 

39 

212.177450 

22.994366 

0.4188 

0.424 ± 0.015 

25.4 

2.28 

l.I 

0.4522 

13.1 

11.9 

40 

213.120494 

22.548543 


0.420 ± 0.014 

38.4 

3.57 

1.3 

0.6800 

19.8 

15.4 

41 

211.935365 

23.227769 

0.4199 

0.416 ± 0.019 

10.7 

0.90 

0.8 

0.6825 

19.9 

24.5 

42 

212.374306 

22.818205 


0.420 ± 0.019 

17.3 

1.50 

1.0 

0.3710 

10.8 

11.2 

43 

212.420269 

23.093511 


0.428 ± 0.018 

14.5 

1.24 

0.9 

0.2209 

6.6 

7.3 

44 

212.593164 

22.699892 


0.430 ± 0.017 

12.1 

1.02 

0.8 

0.3848 

11.5 

13.6 

45 

211.957129 

23.106715 


0.434 ± 0.014 

28.6 

2.59 

1.1 

0.6471 

19.4 

16.9 

46 

212.999005 

22.821469 

0.4358 

0.443 ± 0.019 

12.5 

1.06 

0.9 

0.4052 

12.2 

14.3 

47 

212.607221 

22.800250 


0.440 ± 0.018 

13.0 

1.10 

0.9 

0.2837 

8.6 

10.0 

48 

212.115194 

23.041513 

0.4397 

0.442 ± 0.016 

11.1 

0.94 

0.8 

0.5027 

15.3 

18.7 

49 

212.968630 

22.717594 


0.441 ± 0.021 

11.0 

0.92 

0.8 

0.4604 

14.0 

17.3 

50 

212.980494 

22.615241 


0.441 ± 0.018 

15.6 

1.34 

0.9 

0.5505 

16.7 

18.2 

51 

212.172529 

23.471857 


0.445 ± 0.018 

15.6 

1.35 

0.9 

0.5952 

18.3 

19.8 

52 

212.710891 

23.162554 


0.452 ± 0.017 

11.4 

0.96 

0.8 

0.0951 

3.0 

3.6 

53 

212.336805 

22.647284 


0.454 ± 0.019 

15.5 

1.34 

0.9 

0.5251 

16.4 

18.0 

54 

212.658634 

23.039911 

0.4582 

0.437 ± 0.014 

37.2 

3.44 

1.3 

0.0399 

1.3 

1.0 

55 

212.346270 

23.242090 

0.4585 

0.449 ± 0.016 

15.6 

1.35 

0.9 

0.3310 

10.4 

11.4 

56 

213.213654 

22.988852 

0.4603 

0.477 ± 0.013 

47.9 

4.52 

1.4 

0.5176 

16.4 

11.9 

57 

213.210168 

23.112030 

0.4615 

0.469 ± 0.013 

48.9 

4.62 

1.4 

0.5071 

16.1 

11.7 


(1) Cluster ID. (2) Right ascension (J2000). (3) Declination (J2000). (4) Spectroscopic redshift. (5) Photometric redshift. (6) Richness of galaxies 
having L ^ 0.2L*, corrected for incompleteness (hence non-integer). (7) Inferred mass using Equation |T]; typical scatter of ~ 0.25 in ln(M) 
jRykoff et al.p012j . (8) Inferred virial radii of the cluster using Equation (9) Projected separation to the Q1410 sightline in degrees. (10) Projected 
separation to the Q1410 sightline in Mpc. (11) Projected separation to the Q1410 sightline in units of our J2200 estimation. 


imposed a magnitude limit of r < 17.5 mag to select relatively 
bright QSOs. We gave priority to QSOs 2 > 0.3, ensuring large 
redshift path coverage. We also searched in the Galaxy Evolution 
Explorer (GALEX) (Martin et al.|2005l l database and prioritized 
those QSOs with high FUV fluxes to ensure no higher -2 Lyman 
Limit Systems (LLS) were presentj^enabling a signal-to-noise 

^ We believe that biasing against LLS is unimportant for our present 
study. 


ratio (S/N) S/N ~ 10 spectra to be observed in a relatively short 
exposure time (no larger than 15 orbits). For each of these QSOs 
we counted the number of independent cluster-pairs (defined as 
those which were separated by more than 1000 km s“^ from an¬ 
other, and by more than 5000 km s“^ from the background QSO 
in rest-frame velocity space) at impact parameters Ad ^ 2 Mpc 
from the QSO sightline (see Figure [T] for an illustration). We 
note that in this paper we adopted a larger limit of 3 Mpc as the 
fiducial minimum impact parameter (see Section]^, motivated by 
the results obtained in Section|^ We finally selected the sightline 
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that maximized the number of independent cluster-pair structures, 
for the mininum observing time. Table[T]presents a summary of 
the properties of the targeted QSO, Q1410. 

2.2.2 HST/COS observations and data reduction ofQMlO 

In this paper we present FUV spectroscopic data of Q1410 (see 
Table for a summ ary of its main properties) from HST/COS 
(|Green et al.|2012J taken under program General Observer (GO) 
12958 (PI Tejos). 

The QSO QI4I0 was observed in August 2013 using both 
GI30M and GI60M gratings centred at 1318 and 1611 A re¬ 
spectively, using the four fixed-pattern noise positions (FP-POSs) 
available for each configuration. These settings provided medium- 
resolution (i? = ^ ~ 16000 — 21000 or FWHM ~ 
0.07—0.09A)overthe FUV wavelength range of ~ 1160—1790 
A, but having two ~ 20 A gaps around each central wavelength. 
We chose this approach rather than a continuous wavelength cov¬ 
erage to increase the S/N at spectral regions where we expected 
to find H I absorption line systems associated with inter-cluster 
filaments. Having multiple FP-POS (as in our observations) is 
crucial, however, to minimize the effect of fixed-pattern noise 
present in COS (see |Osterman et al.|201 l[|Green et al.|2012| for 
more details on the technical aspects of COS). 

Data reduction was performed in the same fashion as pre¬ 
sented in |Finn et ^^ ( [2014^ and |Tejos et al.| l |2014) , for the G130M 
and G160M COS gratings. In summary, we used the CALCOS 
v2.I8.5 pipeline with extraction windows of 25 and 20 pixels for 
the G130M and G160M gratings, respectively. We applied a cus¬ 
tomized background estimation smoothing (boxcar) over 1000 
and 500 pixels for the FUVA and FUVB stripes, respectively, 
while masking out and linearly interpolating over regions close 
to geocoronal emission lines and pixels flagged as having bad 
quality. The uncertainty was calculated in the same manner as in 
CALCOS but using our customized background. The co-alignment 
was performed using strong Galactic interstellar medium (ISM) 
features as reference. We finally binned the original spectra hav¬ 
ing dispersions of ~ 0.010 and ~ 0.012 A pixeU^ for the 
G130M and G160M gratings respectively, into a single linear 
wavelength scale of 0.0395 A pixeU^ (roughly corresponding 
to two pixels per resolution element). Due to the difficulties in 
assessing the degree of geocoronal contamination in the final re¬ 
duced Q1410 spectrum, we opted to mask out the spectral regions 
close to rest-frame N I, H I Lya and O I (namely 1300.0—1307.5, 
1198.5 - 1201.0 and 1213.5 - 1217.8 A, respectiv ely). 

Our pseudo-continuunj^fit was modelled as in |Tejos et al. 
P014^ , but also introducing the presence of thre e partial LLS 
breaks at « 1232, 1401 and 1637 A. Figure shows the re¬ 
duced Q1410 spectrum (black line), its corresponding uncertainty 
(green lines) and our adopted pseudo-continuum fit (blue dotted 
line). 

3 CHARACTERIZATION OF LARGE-SCALE 

STRUCTURES AROUND THE Q1410 SIGHTLINE 

3.1 Galaxy clusters 

From the redMaPPer catalog described in Section [ZT] we define 
a subsample of clusters according to the following criteria: 

• the redshift has to lie between 0.1 ^ z ^ 0.5; and. 


® i.e. including intrinsic broad emission lines and the Galaxy’s damped 
Lya system wings. 


• the impact parameter to the Q1410 sightline has to be no 
larger than 20 Mpc. 

There are a total of 57 clusters from the redMaPPer catalog 
satisfying the aforementioned criteria, whose relevant informa¬ 
tion is presented in Table We also show their distribution 
around the Q1410 sightline in Figure|^(coloured circles). 

The redshift range of 0.1 ^ 2 ^ 0.5 was chosen to ensure 
simultaneous coverage of both H I and O VI transitions from our 
COS data, while the impact parameter of 20 Mpc (arbitrary) was 
chosen to cover scales expected to be relevant for inter-cluster 
filaments (e.g. |Colberg et al.|2005[|Gonzalez & Padilla|2010) . 

In Appendix |a| we show how our subsample of clusters 
compares to appropriate control samples drawn from the full 
redMaPPer catalog. We found no statistically significant differ¬ 
ences for the mass (richness) and redshift distributions between 
our subsample and the control samples, implying that no notice¬ 
able bias is present in the subsample close to the Q1410 sightline. 


3.2 Cluster-pairs 

From the subsample of clusters around the Q1410 sightline pre¬ 
sented in Table|^ we define a sample of cluster-pairs according 
to the following criteria: 

• the rest-frame velocity difference between the clusters red- 
shifts has to be < 2000 km s“^; 

• at least one of the two members has to have a spectroscopic 
redshift determination (typically from a BCG), and the other has 
to have a redshift uncertainty no larger than 0.05[^ 

• the transverse separation between the cluster centres has to 
be no larger than 25 Mpc; and, 

• the impact parameter between the inter-cluster axis and the 
Q1410 sightline has to be Ad < 3 Mpc. 

When these criteria are satisfied, we assign the cluster-pair 
redshift to be the average between the two cluster members. 
There are a total of 11 cluster-pairs satisfying these criteria 
around the Q1410 sightline (see the grey dashed lines in Fig¬ 
ure]^, and whose relevant information is presented in Table[^ 
We chose 2000 km s“^ (arbitrary) for the rest-frame ve¬ 
locity difference limit for the clusters in a cluster-pair, in order 
to account for the typical velocity dispersion of galaxy clus¬ 
ters (~ 600 km s“^) and a contribution from a cosmological 
redshift difference. We note however that the majority of the clus¬ 
ters in a given cluster-pair have rest-frame velocity differences 
< 1000 km s“^ and that all of them have < 1400 km s“^ (see 
the fifth column of Table[^. The need for relatively small redshift 
uncertainties for the clusters is necessary to minimize the dilution 
of a real signal due to unconstrained positions for the cluster- 
pairs along the line-of-sight. The 25 Mpc (arbitrary) maximum 
separation between clusters in a cluster-pair was motivated by 
theoretical results from W-body simulations in ACDM universes. 
These studies indicate that at < 25 Mpc there is relatively high 
probability of having coherent filamentary structures between 
galaxy clusters (e.g. |Colberg et al.|2005[ [Gonzalez & Padilla| 
|2010| l. We stress that the majority of cluster-pairs in our sample 
have projected separations ~ 10—15 Mpc (see the fourth column 
of Table[^. The choice for the maximum impact parameter be¬ 
tween the cluster-pair inter-cluster axis and the Q1410 sightline 


® Although a photometric uncertainty of 5z = 0.05 corresponds to 
a very large &v 15 000/(1 -I- z) km s“^, we note that in most of 
our cluster-pair sample both clusters have spectroscopic redshifts (see 
Table[^. 
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Figure 2. Representation of the LSSs within 20Mpc around the Q1410 sightline at redshifts 0.1 ^ 2 ^ 0.5, as traced by the 57 galaxy clusters from 
the redMaPPer catalog (see Section[^and Table[^. Galaxy clusters are represented by coloured circles according to redshift, as given by the colour 
bar scale. Cluster-pairs at impact parameters Ad ^ 3Mpc of the Q1410 sightline are represented by dashed grey lines. The left panel shows to the 
projected in the sky distributions in co-moving Mpc, where the Q1410 sightline is represented by the yellow star at the origin. The right panel shows the 
distribution along the line-of-sight in terms of total transverse separation (in co-moving Mpc) as a function of redshift (note that inter-cluster axes appear 
as hyperbolas in these coordinates). See Section[^for further details. 


Table 3. Clusters-pairs around Q1410. 


Pair ID 

Cluster IDs 

z 

Separation between clusters 

Ad 

Arc 

Both spec-z? 

Grouped ID 




Transverse (Mpc) 

LOS (knis“4) 

(Mpc) 

(Mpc) 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

1 

9,10 

0.1588 

14.7 

419 

0.48 

0.31 

y 

1 

2 

12,13 

0.1760 

14.7 

1378 

1.54 

7.33 

y 

2 

3 

23,25 

0.3439 

12.7 

1189 

1.85 

3.43 

y 

3 

4 

25,28 

0.3488 

18.1 

1021 

2.29 

9.00 

y 

4 

5 

34,35 

D3126 

17.3 

771 

2.76 

4.63 

n 

5 

6 

37,42 

0.4139 

13.2 

1365 

2.59 

2.79 

n 

6 

7 

37,38 

0.4149 

5.3 

435 

1.86 

1.99 

y 

6 

8 

37,39 

0.4163 

16.7 

1043 

1.16 

3.62 

y 

6 

9 

37,41 

0.4169 

23.6 

1284 

0.03 

3.80 

y 

6 

10 

54,55 

0.4584 

11.1 

68 

1.05 

0.69 

y 

7 

11 

54,57 

0.4599 

16.2 

682 

1.24 

0.21 

y 

7 


(1) Cluster-pair ID. (2) IDs of clusters defining the cluster-pair as given in Table[^ (3) Redshift of the cluster-pair. (4) Transverse 
separation between clusters in Mpc. (5) Along the line-of-sight separation between the clusters in rest-frame km s~ ^. (6) Impact 
parameter to the Q1410 sightline. (7) Projected on the sky distance to the c/ose.?r cluster of the pair, along the inter-cluster axis. (8) 
Whether both clusters have spectroscopic redshifts. (9) Grouped ID for independent cluster-pairs. 


of Ad = 3 Mpc was directly motivated by one of our obser¬ 
vational results (see Section|^, and is in good agreement with 
the typical scales for the radii of inter-cluster filaments inferred 
from A'^-body simulations (e.g.|Colberg et al.|2005||Gonzalez &| 
|Padilla|20 1 0[ | Aragon-Calvo et al. |20 1 OlT 

3.3 Independent cluster-pairs 

As expected from the clustering of galaxy clusters, many cluster- 
pairs are grouped together and hence might not be tracing inde¬ 
pendent structures. We therefore have grouped cluster-pairs if 
they are within 1000 km s“^ from one another and we treat them 
as independent. There are a total of 7 independent cluster-pairs; 
a unique identifier is given for each of these in the last column of 
Table We can clearly see these structures in the right panel of 
Figur^at 2 ~ 0.16,0.18, 0.34,0.35,0.37, 0.41 and 0.46. 

Again, this velocity limit of 1000km s“^ is arhitrary and 
chosen to account for the typical velocity dispersion of galaxy 
clusters and a contribution from a cosmological redshift differ¬ 
ence. As reference, if 3000 km s“^ is used instead, then there 
are 6 independent structures rather than 7 (i.e. the structures at 
2 ~ 0.34 — 0.35 are joined together). We note however that in 
our subsequent analysis of associating IGM absorption lines with 


cluster-pairs, we will only use the impact parameter to the closest 
cluster-pair independently of the group it belongs to (see Sec- 
tion|^. Therefore, the velocity limit for grouping cluster-pairs 
is irrelevant for the main results of this paper. Still, this defini¬ 
tion allows us to quantify how many independent cluster-pairs 
the Q1410 sightline is probing, making sure that our results are 
not dominated by a single coherent structure spanning a large 
redshift range. 

3.4 How unusual is the Q1410 field? 

As explained in Section [2.2.1| the field around Q1410 was se¬ 
lected to maximize the presence of cluster-pairs close to the QSO 
sightline. Therefore, it is by no means a randomly selected sight¬ 
line. To quantify how unusual the sightline is, we have performed 
the same search for clusters, cluster-pairs and independent cluster- 
pairs in 1000 randomly selected sightlines having coordinates 
R.A. £ [140, 222] degrees and Dec. £ [4, 56] degrees (i.e. well 
within the SDSS footprint), using the same set of criteria used 
to characterize the Q1410 field (see Section [T7| Section [T2[ and 
Section [T^ . 

In Figure[^we compare our observed (dashed black vertical 
lines) number of redMaPPer clusters (left panel), cluster-pairs 
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Figure 3. Comparison between our observed (dashed black vertical lines) number of redMaPPer clusters (left panel), cluster-pairs (middle panel) and 
independent cluster-pairs (right panel) in the Q1410 sightline satisfying our criteria (see Section [x7] Section [X^ and Section [T^ respectively), and the 
distributions from control samples (solid red histograms; see Section [X4} . The Q1410 sightline is highly exceptional in terms of number of LSS traced by 
galaxy clusters close to it. 


(middle panel) and independent cluster-pairs (right panel), to the 
distributions from our control samples (solid red histograms). 
There are a total of 57 clusters at redshifts 0.1 ^ 2 ^ 0.5 satisfy¬ 
ing the condition of being at impact parameter of < 20 Mpc from 
the Q1410 sightline, whereas the average random expectation is 
32±10 with median of 32. Likewise, the actual number of cluster- 
pairs and independent cluster-pair within our constraints are 11 
and 7 respectively, whereas the average random expectations are 
1.6 ± 2.2 and 1 ± 1, with medians of 1 and 1, respectively. These 
last two distributions are very skewed towards zero. 

Although the number of clusters around Q1410 exceeds 
that from the random expectation at only the 2 — 3a confidence 
level (c.l.), the excesses of total and independent cluster-pairs 
are highly significant (> 5a), making Q1410 a very exceptional 
sightline. We take this fact into account when comparing the 
incidences of absorption line systems close to cluster-pairs and 
the field as estimated from the Q1410 sightline itself (see Sec- 
tion |5.2[ (. 


4 CHARACTERIZATION OF ABSORPTION LINES IN 
THE Q1410 SPECTRUM 

We performed a. full characterization of absorption lines in the 
HST/COS FUV spectrum of Q1410. This approach is more time 
consuming than just restricting ourselves to spectral regions as¬ 
sociated with the redshifts where known structures exist (e.g. 
clusters, cluster-pairs; see Section]^, but is necessary to avoid 
potential biases and systematic effects. In particular, our approach 
allows us to: (i) identify absorption lines independently of the 
presence of known structures; (ii) quantify how the rest of the 
redshift path unassociated with these known structures compares 
to the field expectation in terms of absorption features (see Sec- 
tion |5.2| l; and (iii) assess the extent of contamination by blended 
unassociated lines in a given redshift. This last point is crucial to 
minimize misidentification of lines, but some ambiguous cases 
are unavoidable. In this section we present our methodology for 
the identification and characterization of absorption lines in the 
Q1410 spectrum, and how we handled ambiguity. 


4.1 Absorption line identification 

We searched for individual absorption line component j^in the 
continuum normalized QSO spectrum manually (i.e. eyeballing), 
based on an iterative algorithm described as follows; 

(i) Identify all possible absorption components (H I and met¬ 
als) within ±500 km s“^ from redshift 2 = 0, and assign them 
to the ‘reliable’ category (label ‘a’; see Section [43) l. 

(ii) Identify all possible absorption components (H I and met¬ 
als) within ±500 km s”'^ from redshift z = 2qso, and assign 
them to the ‘reliable’ category. 

(iii) Identify H I absorption components, showing at least two 
transitions (e.g. Lya and Ly/3 or Ly/3 and Ly 7 , and so on; i.e. 
strong H l|^ starting at 2 = 2qso until 2 = 0, and assign them 
to the ‘reliable’ category. This identification includes the whole 
Lyman series covered by the spectrum in a given component. 

(iv) Identify all possible metal absorption components within 
±200 km s“^ from each H I redshift found in the previous step, 
and assign them to the ‘reliable’ category. When the wavelength 
coverage allows the detection of multiple transitions of a single 
ion, we require the relative positions of these to coincide; in the 
case of multiple adjacent components blending with each other, 
we require them to have similar kinematic structure across the 
multiple transitions of the same ion. 

(v) Identify high-ionization transitions (namely: Ne VIII, O VI, 
N V, C IV, Si IV) showing in at least two transitions, indepen¬ 
dently of the presence of H I, starting at 2 = 2qso until 2 = 0, 
and assign them to the ‘reliable’ category. 

(vi) Identify low-ionization transitions (namely: C II, C III, 
N II, N III, O I, O II, Si II, Si III, Fe II, Fe III and Al II), showing 
at least two transitions, independently of the presence of H I, 
starting at 2 = 2qso until 2 = 0, and assign them to the ‘reliable’ 
category]^ 


In this paper a ‘component’ is defined as an individual absorption 
line in a given ion; and a ‘system’ is defined as the union of components 
lying within a given velocity window (usually arbitrary) from a given 
redshift. 

Note that this condition allow us to identify strong H I components at 
redshifts greater than z > 0.477 by means of higher order Lyman series 
transitions. 

No low-ionization transition was found without having associated 
H I, so this step was redundant. 
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Table 4. Absorption lines in the Q1410 sightline. 


Component ID 

Ion 

Obs. wavelength 

(A) 

z 

log(N’/cm 

b 

(km s~^) 

Wr 

(A) 

(S/iV>„s 

Label 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


1 

C IV 

1547.8 

-0.00024 ± 0.00013 

13.98 ± 0.82 

18 ±35 

0.205 ± 0.383 

14 ± 1 

a 

2 

Si IV 

1393.4 

-0.00024 ± 0.00003 

13.38 ± 0.82 

15 ± 22 

0.134 ± 0.231 

12 ± 1 

a 

3 

A1II 

1670.5 

-0.00018 ± 0.00100 

14.13 ± 40.72 

11 ± 131 

0.268 ± 2.683 

11 ± 1 

a 

4 

N V 

1238.6 

-0.00018 ± 0.00011 

13.55 ± 0.53 

23 ±55 

0.064 ±0.105 

15 ± 1 

a 

5 

Si IV 

1393.6 

-0.00014 lii 0.00008 

13.73 ± 0.38 

35 ± 14 

0.304 ± 0.203 

12 ± 1 

a 

6 

Si III 

1206.4 

-0.00012 ± 0.00003 

16.24 ± 1.53 

16 ± 10 

0.534 ± 1.077 

14 ± 1 

a 

7 

C IV 

1548.0 

-0.00011 ± 0.00011 

14.09 ± 0.66 

18 ± 29 

0.227 ± 0.340 

15 ± 1 

a 

8 

Si II 

1260.3 

-0.00010 ± 0.00001 

16.58 ± 0.08 

17 ± 1 

0.700 ± 0.057 

16 ± 1 

a 

9 

C II 

1334.4 

-0.00010 ± 0.00003 

18.32 ± 1.65 

23 ±9 

0.824 ± 1.500 

11 ± 1 

a 

10 

A111 

1670.7 

-0.00006 ± 0.00107 

14.07 ± 66.61 

9 ± 204 

0.222 ± 2.217 

11 ± 1 

a 

11 

C II* 

1335.6 

-0.00006 ± 0.00007 

14.11 ± 0.38 

24 ±32 

0.164 ±0.188 

12 ± 1 

a 

12 

Si III 

1206.7 

0.00019 ± 0.00004 

13.28 ± 3.99 

11 ±41 

0.121 ±0.894 

14 ± 1 

a 

13 

Si III 

1207.2 

0.00056 ± 0.00016 

12.37 ± 0.92 

17 ± 88 

0.043 ±0.169 

13 ± 1 

a 

14 

Hi 

1222.7 

0.00579 ± 0.00003 

12.86 ± 0.17 

19 ± 17 

0.036 ±0.019 

13 ± 1 

c 

15 

Hi 

1224.7 

0.00744 ± 0.00002 

12.90 ± 0.12 

12 ± 11 

0.037 ± 0.020 

14 ± 1 

c 

16 

Hi 

1225.0 

0.00771 ± 0.00002 

12.92 ± 0.30 

7 ± 14 

0.034 ± 0.032 

14 ± 1 

c 

17 

Hi 

1250.5 

0.02865 ± 0.00001 

13.49 ± 0.06 

18 ±5 

0.115 ± 0.022 

16 ± 1 

b 

18 

Hi 

1250.9 

0.02894 ± 0.00001 

13.57 ± 0.05 

32 ±7 

0.155 ±0.023 

16 ± 1 

b 

19 

Hi 

1259.4 

0.03594 ± 0.00001 

13.74 ± 0.05 

27 ± 4 

0.197 ±0.025 

16 ± 1 

b 

20 

Hi 

1268.9 

0.04375 ± 0.00001 

13.49 ± 9.56 

4 ± 26 

0.044 ± 0.441 

15 ± 1 

c 

21 

Hi 

1269.4 

0.04418 lii 0.00004 

12.93 ± 0.16 

30 ± 19 

0.043 ±0.017 

15 ± 1 

c 

22 

Hi 

1292.0 

0.06277 ± 0.00001 

13.32 ± 0.08 

21 ±7 

0.089 ±0.019 

14 ± 1 

b 

23 

Hi 

1295.9 

0.06599 ± 0.00004 

12.86 ± 0.17 

21 ± 16 

0.036 ±0.017 

14 ± 1 

c 

24 

Hi 

1298.3 

0.06800 ± 0.00077 

14.32 ± 6.43 

25 ± 76 

0.300 ± 3.003 

14 ± 1 

c 

25 

Hi 

1324.8 

0.08981 ± 0.00002 

13.29 ± 0.10 

22 ± 10 

0.087 ± 0.025 

10 ± 1 

b 

26 

Hi 

1372.7 

0.12920 ± 0.00002 

13.55 ± 0.06 

44 ± 8 

0.159 ±0.023 

12 ± 1 

b 

27 

Hi 

1410.9 

0.16058 ±0.00005 

13.40 ± 0.11 

59 ± 22 

0.123 ±0.033 

14 ± 1 

b 

28 

Hi 

1416.4 

0.16515 ± 0.00006 

12.80 ± 0.28 

23 ± 26 

0.032 ± 0.026 

14 ± 1 

c 

29 

Hi 

1419.3 

0.16748 ± 0.00005 

12.76 ± 0.17 

30 ± 19 

0.030 ±0.013 

16 ± 2 

c 

30 

Hi 

1424.7 

0.17192 ± 0.00004 

12.88 ± 0.13 

34 ± 16 

0.039 ±0.013 

19 ±2 

b 

31 

Hi 

1429.4 

0.17580 ± 0.00004 

12.86 ± 0.15 

31 ± 17 

0.037 ±0.014 

22 ± 2 

c 

32 

Hi 

1433.2 

0.17894 ± 0.00002 

13.07 ± 0.09 

23 ± 8 

0.057 ±0.013 

21 ± 4 

b 

33 

Hi 

1441.1 

0.18543 ± 0.00002 

13.00 ± 0.07 

22 ± 6 

0.048 ± 0.009 

24 ± 1 

b 

34 

Hi 

1444.1 

0.18794 ± 0.00007 

13.41 ± 0.08 

112 ± 28 

0.133 ±0.026 

23 ±2 

b 

35 

Hi 

1449.9 

0.19267 ± 0.00002 

12.91 ± 0.09 

24 ±9 

0.041 ± 0.009 

23 ± 2 

b 

36 

C 11 

1613.5 

0.20904 ± 0.00001 

14.53 ± 1.73 

6±6 

0.096 ±0.135 

8 ± 2 

a 

37 

Si II 

1523.9 

0.20905 ± 0.00002 

12.54 ± 0.10 

11 ± 8 

0.046 ±0.016 

16 ± 1 

a 

38 

Hi 

1469.8 

0.20908 ± 0.00001 

16.15 ± 1.03 

11 ± 3 

0.257 ±0.155 

18 ± 1 

a 

39 

Si III 

1458.8 

0.20909 ± 0.00001 

12.63 ± 0.06 

13 ± 4 

0.066 ± 0.011 

20 ± 1 

a 

40 

C III 

1181.4 

0.20915 ± 0.00001 

14.60 ± 13.57 

8 ±35 

0.113 ± 1.134 

12 ± 1 

a 

41 

N III 

1196.8 

0.20916 ± 0.00002 

13.84 ± 0.10 

19 ± 10 

0.060 ±0.017 

13 ± 1 

a 

42 

Hi 

1470.0 

0.20920 ± 0.00030 

13.37 ± 2.72 

14 ±47 

0.088 ±0.571 

18 ± 1 

a 

43 

Si II 

1524.5 

0.20955 ± 0.00001 

13.05 ± 0.15 

13 ±5 

0.109 ± 0.043 

15 ± 2 

a 

44 

C 11 

1614.2 

0.20955 ± 0.00001 

14.96 ± 1.42 

10 ±7 

0.174 ± 0.184 

8 ± 2 

a 

45 

Si IV 

1685.8 

0.20956 ± 0.00001 

13.67 ± 0.08 

14 ±3 

0.170 ±0.033 

10 ± 1 

a 

46 

Hi 

1470.5 

0.20958 ± 0.00015 

16.29 ± 3.61 

14 ± 7 

0.330 ± 3.272 

18 ± 1 

a 

47 

N III 

1197.3 

0.20962 ± 0.00001 

14.94 ± 1.48 

14 ± 12 

0.151 ± 0.192 

13 ± 1 

a 

48 

Si III 

1459.4 

0.20962 ± 0.00001 

13.56 ± 0.07 

23 ±3 

0.256 ± 0.032 

20 ± 1 

a 

49 

C III 

1181.8 

0.20964 ± 0.00236 

15.03 ±59.37 

13 ±98 

0.188 ± 1.884 

12 ± 1 

a 

50 

N V 

1498.5 

0.20964 ± 0.00003 

13.57 ± 0.09 

31 ± 10 

0.071 ± 0.015 

17 ± 1 

a 

51 

O VI 

1248.3 

0.20967 ± 0.00005 

14.31 ± 0.21 

26 ± 10 

0.164 ± 0.073 

16 ± 1 

a 

52 

Si II 

1524.8 

0.20973 ± 0.00005 

13.06 ± 0.16 

34 ± 13 

0.150 ±0.058 

15 ±2 

a 

53 

C II 

1614.5 

0.20976 ± 0.00001 

14.48 ± 0.04 

20 ± 3 

0.233 ± 0.029 

8 ± 2 

a 

54 

Hi 

1470.7 

0.20980 ± 0.00024 

15.93 ± 2.00 

27 ± 15 

0.534 ± 0.567 

18 ± 1 

a 

55 

Si IV 

1686.2 

0.20981 ± 0.00002 

13.27 ± 0.07 

26 ±7 

0.128 ± 0.025 

10 ± 1 

a 

56 

Si III 

1459.7 

0.20986 ± 0.00001 

13.17 ± 0.05 

21 ± 4 

0.178 ±0.025 

20 ± 1 

a 

57 

C III 

1182.1 

0.20986 ± 0.00021 

17.15 ± 0.78 

14 ± 18 

0.488 ±0.515 

12 ± 1 

a 

58 

O VI 

1248.5 

0.20987 ± 0.00006 

14.16 ± 0.31 

26 ± 16 

0.129 ± 0.095 

16 ± 1 

a 

59 

N III 

1197.5 

0.20987 ± 0.00003 

14.09 ±0.11 

25 ± 12 

0.100 ± 0.031 

13 ± 1 

a 

60 

Hi 

1471.3 

0.21028 ± 0.00001 

13.70 ± 0.05 

19 ±3 

0.157 ±0.020 

18 ± 1 

a 

61 

O VI 

1248.9 

0.21028 ± 0.00001 

14.50 ± 0.03 

50 ±6 

0.270 ± 0.022 

16 ± 1 

a 

62 

N V 

1499.3 

0.21029 ± 0.00006 

13.53 ± 0.15 

50 ± 24 

0.068 ± 0.024 

16 ± 1 

a 

63 

C III 

1182.5 

0.21033 ±0.00002 

13.09 ± 0.18 

11 ± 11 

0.055 ± 0.040 

12 ± 1 

a 

64 

Hi 

1478.4 

0.21615 ±0.00002 

13.30 ± 0.05 

42 ±7 

0.096 ±0.012 

17 ± 1 

b 

65 

Hi 

1481.9 

0.21900 ± 0.00002 

13.10 ± 0.07 

36 ±9 

0.063 ±0.011 

17 ± 1 

b 

66 

Hi 

1486.8 

0.22300 ± 0.00002 

13.32 ± 0.05 

52 ± 8 

0.104 ± 0.012 

17 ± 1 

b 


(1) Absorption component ID. (2) Ion (see Section |4T| for details on the line identification process). (3) Observed wavelength of the sti'ongest transition 
of the ion in the HST/COS spectrum. (4) Redshift from the Voigt profile fitting (see Sectionj^j. (5) Column density from the Voigt profile fitting (see 
Section |4^ . (6) Doppler parameter from the Voigt profile fitting (see Section [42) . (7) Inferred rest-frame equivalent width from fitted values (note 
that uncertainties are greatly overestimated for saturated lines or unconstrained fits; see Sectionj^j. (8) Averaged local S/N (see Sectionj^J. (9) Line 
reliability flag (“a’ secure, ‘b’ possible and ‘c’ uncertain; see Section[T4j. 
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Table 4 - continued 


Component ID 

Ion 

Obs. wavelength 

(A) 

z 

log(iV/cm 

h 

(km s“^) 

Wr 

(A) 

(S/AT)™ 

Label 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


67 

Hi 

1497.1 

0.23150 ±0.00001 

13.45 ± 0.04 

32 ±4 

0.125 ± 0.012 

17 ± 1 

b 

68 

O VI 

1274.5 

0.23508 ±0.00001 

14.08 ± 0.05 

24 ±4 

0.111 ± 0.014 

15 ± 1 

a 

69 

Hi 

1501.5 

0,23511 ± 0.00001 

14.48 ± 0.04 

32 ± 1 

0.395 ± 0.021 

17 ± 1 

a 

70 

Hi 

1509.0 

0.24126 ±0.00003 

13.29 ± 0.06 

54 ± 10 

0.097 ± 0.014 

16 ± 1 

b 

71 

Hi 

1538.8 

0.26584 ± 0.00005 

14.91 ± 0.41 

28 ± 4 

0,419 ± 0.111 

15 ± 1 

a 

72 

C III 

1236.8 

0,26587 ±0.00001 

13.70 ±40.19 

3 ±66 

0.041 ± 0.406 

15 ± 1 

a 

73 

Hi 

1539.0 

0,26593 ± 0.00001 

15.88 ± 0.06 

16±2 

0,328 ± 0.041 

15 ± 1 

a 

74 

Hi 

1541.8 

0.26827 ±0.00010 

12.78 ± 0.22 

50 ± 37 

0.032 ± 0.017 

15 ± 1 

c 

75 

Hi 

1545.5 

0.27132 ±0.00004 

12.98 ± 0.12 

36 ± 16 

0.049 ± 0.015 

14 ± 2 

b 

76 

Hi 

1550.9 

0,27574 ± 0.00004 

13.87 ± 0.32 

21 ± 17 

0.200 ± 0.166 

14 ± 1 

c 

77 

Hi 

1566.7 

0,28874 ± 0.00007 

13.06 ± 0.15 

49 ±25 

0.059 ± 0.022 

13 ± 1 

c 

78 

Hi 

1573.8 

0.29459 ±0.00017 

13.58 ± 0.13 

157 ±57 

0.193 ± 0.058 

14 ± 1 

b 

79 

Hi 

1576.7 

0.29696 ±0.00010 

12.75 ± 0.23 

45 ±37 

0.030 ± 0.017 

13 ± 1 

c 

80 

Hi 

1578.8 

0,29874 ± 0.00001 

13.68 ± 0.03 

36 ±4 

0.191 ± 0.015 

13 ± 1 

b 

81 

Hi 

1582.7 

0,30188 ±0.00002 

13.25 ± 0.07 

32 ±8 

0.084 ± 0.014 

13 ± 1 

b 

82 

Hi 

1584.2 

0.30314 ±0.00003 

13.32 ± 0.07 

47 ± 11 

0.102 ± 0.017 

13 ± 1 

b 

83 

Hi 

1613.9 

0.32756 ± 0.00003 

13.24 ±0.11 

24 ± 11 

0.079 ± 0.024 

8 ± 2 

b 

84 

Hi 

1617.4 

0.33044 ± 0.00004 

13.38 ± 0.09 

45 ± 13 

0.113 ± 0.024 

10 ± 1 

b 

85 

Hi 

1631.6 

0.34217 ±0.00006 

13.75 ± 0.05 

153 ± 19 

0,276 ± 0.030 

10 ± 1 

b 

86 

O VI 

1392.6 

0.34954 ± 0.00002 

13.77 ± 0.09 

21 ± 8 

0.061 ± 0.015 

11 ± 1 

a 

87 

Si III 

1628.6 

0.34986 ± 0.00001 

13.50 ± 0.03 

45 ± 4 

0,378 ± 0.032 

10 ± 1 

a 

88 

Hi 

1641.0 

0.34986 ± 0.00001 

15.88 ± 0.03 

17 ± 1 

0.344 ± 0.015 

11 ± 1 

a 

89 

O VI 

1393.0 

0.34989 ±0.00001 

14.07 ± 0.06 

19 ± 5 

0.101 ± 0.018 

11 ± 1 

a 

90 

C II 

1801.5 

0.34989 ±0.00001 

13.31 ± 0.36 

5±5 

0.030 ± 0.028 

11 ± 1 

a 

91 

N III 

1336.2 

0.34998 ± 0.00006 

14.04 ±0.11 

50 ±20 

0.102 ± 0.028 

12 ± 1 

a 

92 

Hi 

1641.4 

0.35020 ± 0.00004 

14.29 ± 0.09 

97 ± 10 

0.701 ± 0.116 

11 ± 1 

b 

93 

O VI 

1393.4 

0.35029 ±0.00001 

14.48 ± 0.05 

35 ±5 

0.233 ± 0.029 

12 ± 1 

a 

94 

Hi 

1641.6 

0.35035 ±0.00001 

14.57 ± 0.06 

26 ± 3 

0.346 ± 0.047 

11 ± 1 

a 

95 

Hi 

1642.4 

0.35106 ± 0.00001 

15.43 ± 0.03 

25 ± 1 

0.440 ± 0.015 

11 ± 1 

a 

96 

Hi 

1652.3 

0,35918 ± 0.00002 

13.61 ± 0.05 

31 ±5 

0.163 ± 0.022 

11 ± 1 

b 

97 

Hi 

1658.1 

0.36397 ± 0.00003 

13.52 ± 0.06 

58 ± 10 

0.156 ± 0.021 

11 ± 1 

b 

98 

Hi 

1664.1 

0.36886 ± 0.00005 

13.25 ±0.11 

50 ± 18 

0.089 ± 0.023 

11 ± 1 

b 

99 

Hi 

1666.8 

0.37106 ± 0.00001 

14.06 ± 0.03 

36 ±2 

0.329 ± 0.022 

11 ± 1 

a 

100 

Hi 

1667.6 

0.37176 ± 0.00003 

13.24 ± 0.08 

29 ±9 

0.081 ± 0.016 

11 ± 1 

a 

101 

Hi 

1673.8 

0.37686 ±0.00001 

13.59 ± 0.04 

31 ± 4 

0.159 ± 0.018 

11 ± 1 

a 

102 

O VI 

1426.5 

0.38238 ± 0.00002 

13.43 ± 0.14 

7 ± 12 

0.026 ± 0.016 

21 ± 2 

c 

103 

Hi 

1680.6 

0.38242 ± 0.00005 

12.91 ± 0.16 

27 ± 17 

0.041 ± 0.017 

10 ± 1 

c 

104 

Hi 

1683.8 

0.38508 ±0.00001 

14.25 ± 0.03 

32 ±2 

0.345 ± 0.019 

11 ± 1 

a 

105 

Hi 

1691.0 

0.39097 ± 0.00007 

13.19 ± 0.14 

53 ±25 

0.078 ± 0.027 

10 ± 1 

c 

106 

Hi 

1704.6 

0.40219 ± 0.00012 

13.28 ± 0.20 

79 ±49 

0.097 ± 0.048 

10 ± 1 

c 

107 

Hi 

1704.7 

0.40224 ± 0.00003 

13.07 ± 0.30 

15 ± 14 

0.053 ± 0.046 

10 ± 1 

c 

108 

Hi 

1716.3 

0,41182 ±0.00005 

13.47 ± 0.09 

62 ± 18 

0.142 ± 0.031 

9± 1 

b 

109 

Hi 

1723.3 

0.41758 ±0.00001 

14.38 ± 0.03 

19 ± 1 

0.254 ± 0.014 

9± 1 

a 

110 

Hi 

1726.5 

0.42022 ± 0.00006 

13.25 ± 0.12 

56 ± 20 

0.090 ± 0.025 

10 ± 1 

b 

111 

Hi 

1737.3 

0.42911 ± 0.00008 

13.08 ± 0.15 

55 ±27 

0.062 ± 0.023 

9± 1 

c 

112 

Hi 

1740.4 

0.43167 ±0.00005 

13.73 ± 0.83 

23 ± 14 

0.179 ± 0.213 

10 ± 1 

a 

113 

Hi 

1740.6 

0.43183 ±0.00083 

13.21 ± 2.74 

37 ± 102 

0.079 ± 0.793 

9± 1 

c 

114 

Hi 

1768.4 

0.45466 ± 0.00006 

13.46 ± 0.08 

81 ± 18 

0.143 ± 0.026 

9± 1 

b 

115 

Hi 

1775.9 

0.46085 ± 0.00035 

13.55 ± 0.25 

192 ± 141 

0.185 ± 0.115 

9± 1 

c 

116 

Hi 

1777.7 

0,46235 ± 0.00002 

13.77 ± 0.07 

26 ± 5 

0.199 ± 0.037 

9± 1 

a 

117 

On 

1281.0 

0.53509 ±0.00001 

13.85 ± 0.25 

10 ± 10 

0.041 ± 0.033 

15 ± 1 

a 

118 

Hi 

1574.6 

0.53510 ± 0.00001 

16.52 ± 0.02 

16± 1 

0,274 ± 0.010 

14 ± 1 

a 

119 

C II 

2048.6 

0.53510 ±0.00001 

13.30 ± 0.56 

10 ± 10 

0.034 ±0.051 

14 ± 1 

a 

120 

Oiii 

1078.2 

0.53512 ±0.00001 

14.51 ± 0.25 

13 ± 5 

0.078 ± 0.036 

14 ± 1 

a 

121 

C III 

1499.9 

0.53520 ±0.00001 

13.98 ± 0.25 

19 ±4 

0.177 ± 0.054 

16 ± 1 

a 

122 

Oiv 

1209.3 

0.53524 ± 0.00003 

14.22 ±0.11 

23 ± 10 

0.076 ± 0.023 

13 ± 1 

a 

123 

Hi 

1574.8 

0.53531 ±0.00003 

14.99 ± 0.15 

17 ± 4 

0.180 ± 0.050 

14 ± 1 

a 

124 

Hi 

1717.1 

0.67402 ±0.00001 

14.80 ± 0.04 

31 ± 3 

0.238 ± 0.024 

9± 1 

a 

125 

C III 

1752.5 

0.79372 ± 0.00003 

14.84 ± 7.41 

4±9 

0.065 ± 0.653 

9± 1 

a 

126 

Oiv 

1412.9 

0.79373 ± 0.00002 

14.68 ± 0.09 

20 ±5 

0.128 ± 0.032 

14 ± 1 

a 

127 

Oiii 

1259.8 

0.79376 ±0.00001 

14.41 ± 14.69 

3 ±22 

0.023 ± 0.234 

16 ± 1 

a 

128 

Hi 

1839.9 

0.79378 ± 0.00005 

15.55 ± 0.10 

37 ±6 

0.426 ± 0.077 

16 ± 1 

a 

129 

C III 

1752.7 

0.79396 ± 0.00002 

15.14 ± 8.93 

10 ± 22 

0.151 ± 1.514 

9± 1 

a 

130 

Oiv 

1413.2 

0.79401 ± 0.00002 

15.51 ± 4.99 

9± 15 

0.097 ± 0.831 

14 ± 1 

a 

131 

Oiii 

1260.0 

0.79404 ±0.00001 

14.59 ± 0.20 

12 ± 3 

0.077 ± 0.024 

16 ± 1 

a 

132 

Hi 

1840.2 

0.79404 ± 0.00005 

14.91 ± 0.42 

18 ± 10 

0.176 ± 0.115 

16 ± 1 

a 


(1) Absorption component ID. (2) Ion (see Sectionj^^for details on the line identification process). (3) Observed wavelength of the strongest transition 
of the ion in the HST/COS spectrum. (4) Redshift from the Voigt profile fitting (see Section [T^ . (5) Column density from the Voigt profile fitting (see 
Sectionj^j. (6) Doppler parameter from the Voigt profile fitting (see Sectionj^J. (7) Inferred rest-frame equivalent width from fitted values (note that 
uncertainties are greatly overestimated for saturated lines or unconstrained fits; see Section |43] . (8) Averaged local S/N (see Section|^j. (9) Line reliability 
flag (‘a’ secure, ‘b’ possible and ‘c’ uncertain; see Sectionj^J. 

^: Could be a very broad H I Ly/3 at redshift z = 0.53531 instead, but we cannot confinn it with our current data. 
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(vii) Assume all the unidentified absorption features to be 
H I Lya and repeat step (iv). If metals satisfying the criteria 
in step (iv) exist, assign the component to the ‘reliable’ cate¬ 
gory; otherwise assign the component to the ‘possible’ category 
(label ‘b’; see Section [431 l. 

(viii) For complex blended systems we allowed for the pres¬ 
ence of extra heavily blended (hidden) ions, preferentially 
H I Lya unless a metal ion showing at least one unblended tran¬ 
sition exist, and update the identification accordingly. In cases 
where the metal ion shows at least two unblended transitions, 
assign them to the ‘reliable’ category. In the rest of the cases 
(including H I Lya only) assign them to the ‘possible’ category. 

We note that we will then degrade some of the components in the 
‘possible’ category to the ‘uncertain’ category (‘c’), based on an 
equivalent width significance criterium in Section [431 

For all the identified components, we set initial guesses for 
their redshifts, column densities (N), and Doppler parameters (6), 
which are used as the inputs of our automatic Voigt profile fitting 
process described in the Section [4^ We based these guesses 
on the intensities and widths of the spectral features, keeping 
the number of individual components to the minimum: we only 
added a component when there is a clear presence of multiple 
adjacent local minima or asymmetries. In the case of symmetric 
and intense H I Lya absorption lines showing no corresponding 
H I Ly/3 absorption (when the spectral coverage and the S/N 
would have allowed it), this last condition would require the 
components to have relatively large &-values (typically > 40 — 
50km s“^). We warn the reader that this is a potential source 
of bias, especially for the broad Lya systems ()S 50 km s“^) 
expected to trace portions of the WHIM (but see Section [73^ . 

4.2 Voigt profile fitting 

We fit Voigt profiles to the identified absorption line components 
using VPFnl^ We accounted for the non-Gaussian COS line 
spread function (LSF), by interpolating between the closest COS 
LSF tabulated values provided by the Space Telescope Science 
Institute (STSclf^at a given wavelength. We used the guesses 
provided by the absorption line search (see Section|4^ as the 
initial input of VPFIT, and modified them when needed to reach 
converged solutions with low reduced 

When dealing with H I absorption lines, we used at least 
two spectral regions associated with their Lyman series transi¬ 
tions when the spectral coverage allowed it. This means that for 
those showing only the Lya transition, we also included their 
associated Ly/3 regions (even though they do not show evident 
absorption) when available. This last condition provides reliable 
upper limits to the column density of these components. For 
strong H I components, we used regions associated with as many 
Lyman series transitions as possible, but excluding those heavily 
blended or in spectral regions of poor S/N (< 1 per pixel). For 
metal transitions we used all spectral regions available. 

We fitted absorption line systems starting from z = zqso 
until 2 = 0. When a given system at redshift 0 < 2 2qso 
showed strong blends from lower redshifts, we fitted them all si¬ 
multaneously in a given VPFIT iteration (i.e. including all spectral 
regions associated with them). When a given system at redshift 
0 < 2 < 2qso showed weak blending from higher redshifts, we 

http://www.ast.cam.ac,uk/~rfc/vpfit.html 

http://www.Stsci.edu/hst/cos/performance/ 
spectral_resolution 

Our final reduced have in average (and median) values of 1.1. See 
also residuals in FigurelEll 


allowed VPFIT to modify the ‘effective’ continuum by adding 
the previously found absorption line solutions to it. This last con¬ 
dition accounts for the blending of weak lines (especially from 
higher order Lyman series)—whose solutions are already well 
determined—in a more efficient manner than fitting all regions 
involved simultaneously. In the whole process, we allowed VPFIT 
to add lines automatically when the and the Kolmogorov- 
Smirnov test (K-S) test probabilities were below 0.01 (see VPFIT 
manual for details). 

Table shows our final list of identified absorption line 
components, and their corresponding fits. Unique component 
IDs are given in the first column to components for each ion 
(second column). The observed wavelength associated with the 
strongest transition of an ion is shown in the third column (but 
note that some ions can show up in multiple wavelengths when 
having multiple transitions). The fitted redshifts, column densities 
and Doppler parameters are given in the fourth, fifth and sixth 
columns respectively. Our final reduced have an average (and 
median) of 1.1. In Figure [ET] we show how these fits (red line) 
compare to the observed spectrum (black line), by means of 
their residuals (grey dots) defined as the difference between the 
two (i.e. in the same units as the spectrum). We see how these 
residuals are mostly distributed within the spectrum uncertainty 
level (green line). 

4.3 Rest-frame equivalent widths 

For each component we estimate the rest-frame equivalent width 
of the strongest transition, Wr, using the approximation given 
by |Drain^P011| see his equation 9.27), based on their fitted 
N and b values. The resulting values are given in the seventh 
column of Table We chose this approach in order to avoid 
complications when dealing with blended components. We em¬ 
phasize that passing from Wi —>■ {N, b) is not always robust 
when on of the flat part of the curve-of-growth, but passing from 
{N, b) —/ Wi is robust. We compared the results in Wi from our 
adopted approach and that from a direct pixel integration, in a 
subsample of 10 unblended and unsaturated lines, and obtained 
consistent results within the uncertainties. 

The rest-frame equivalent width uncertainty, SWr, was 
estimated as follows. We first calculated the maximum/mini¬ 
mum equivalent width, stpj consistent within Icr 

from the N and b fitted values, i.e. using the aforementioned 
approximation for (N ± 5N, b ± 5b) —/ where 

5N and Sb are the column density and Doppler parameter 
uncertainties, respectively, as given by VPFIT. We then took 
SWi = |(lUr“ - Wr’^). In catastrophic cases where the 
fits are unconstrained (i.e. SWr Wi), we arbitrarily imposed 
5Wi = 10 X Wi, ensuring a very low significance level. 

By using the actual fitted parameters and their correspond¬ 
ing errors, our 5Wi uncertainty estimation takes into account the 
non-Gaussian shape of the COS LSF (particularly important for 
broad absorption lines). For saturated lines, our method will give 
unrealistically large SWi uncertainties due to a poor constraint 
in N, which is a conservative choice. 

4.4 Local signal-to-noise ratio (S/N) estimation 

For each component we estimated the average local spectral 
signal-to-noise ratio (S/N) per pixel, (5'/A)pixei, over the 50 
closest pixels around its strongest transition, without considering 
those with flux values below 90% of the continuum. We then 
estimated the local signal-to-noise ratio (S/N) per resolution 
element {S/W)res as x/2(<S'/A^)pixei. The resulting values are 
given in the eight column of Table|^ 
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Figure 4. Distribution of Doppler parameters (b) as a function of column 
density (N) for H 1 (black circles) and O VI (green squares) absorption 
lines found in the HST/COS FUV spectrum between 0.1 ^ ^5 0.5 and 

excluding those in category ‘c’ (uncertain; see Section [43} . White stars 
mark absorbers associated to cluster-pairs (see SectionBjTThe horizontal 
dashed grey line corresponds to a Doppler parameter of fe = 50 km s”'^, 
our adopted limit to split H I lines into broad and narrow. 


4.5 Absorption line reliability 

To deal with ambiguity and significance of the absorption lines 
we have introduced a reliability flag scheme as follows: 

• Reliable (‘a Absorption line components showing at least 
two transitions or showing up in at least two ions, independently 
of the significance of its corresponding M4. 

• Probable (‘b’): Absorption line components showing 
only one transition, showing up in only one ion, and having 
w.jm, ^ 3. 

• Uncertain (‘c’): Absorption line components showing 
only one transition, showing up in only one ion, and having 
Wi/5Wr < 3. Components in this category will be excluded 
from the main scientific analysis presented in this paper. 

This reliability scheme applied to our absorption line list 
is shown in the ninth column of Table |4] We also show these 
flags together with the ion component ID given (first column of 
Tablel^l in Figure [ eT| as vertical labels. 

We note that previous studies on broad Lya (BLA) lines 
have suggested one to quantifying the completeness level of these 
broad lines by means of {S/N)res x N/b ^Richter et al.|200^ . 
The motivation of this criterion is that the BLA is sensitive to 
both S/N and the optical depth at the line centres, tq oc N/b. 
Therefore, it is not appropriate to use the commonly adopted 
formalism based on a minimum equivalent width threshold for 
unresolved lines (these broad lines are usually resolved). In Ap¬ 
pendix]^ we compare the proposed approach by [Richter et al.| 
^2006| l to ours, and show that imposing a minimum Wr/5Wr 
value (as defined here) is roughly consistent with imposing a 
minimum {S/N)reB x N/b value for broad linej^ but is more 
conservative when applied to narrow lines. Moreover, our ap¬ 
proach has the advantage of being straightforwardly applicable 


We note that jPanforth et al^j2010) reached a similar conclusion 
although using the commonly adopted formalism based on a minimum 
equivalent width threshold for unresolved lines. 


to any absorption line irrespective of its Doppler parameter and 
ionic transition, hence more appropriate for an homogeneus anal¬ 
ysis. 

Figurej^shows the distribution Doppler parameters 6 as a 
function of column densities N, for our sample of H I (black 
circles) and O VI (green squares) between 0.1 ^ 2 : Si 0.5 and 
excluding those in the ‘c’ category (i.e. uncertain). White stars 
mark components that lie within Av ^ 1000 km s“^ and within 
impact parameters of Ad ^ 3 Mpc from cluster-pairs, our fidu¬ 
cial values for associating absorption lines with cluster-pairs (see 
Section]^. 


5 REDSHIFT NUMBER DENSITY OF ABSORPTION 
LINE SYSTEMS AROUND CLUSTER-PAIRS 

After having characterized LSS traced by galaxy cluster-pairs 
around the Q1410 sightline (Section]^ and intervening absorp¬ 
tion lines (Section]^, we can now provide a cross-match between 
the two. Because the completeness level of the redMaPPer clus¬ 
ters with richness < 20 (< 40) at 2 ~ 0.1 — 0.4 (2 ~ 0.4 — 0.5) 
is lower than ~ 50% (see top panel of fig. 22 of |Rykoff et al.j 
|2014| l, in this paper we will only match absorption lines close 
to the position of known cluster-pairs rather than the other way 
around (or both). We also note that the purity of redMaPPer clus¬ 
ters is fairly constant with richness and redshift (some trends are 
present though), but with values above 95% in all the cases (see 
bottom panel of fig. 22 of jRykoff et al.|20T4| >; still, the presence 
of fake clusters will only dilute any real signal when associating 
absorption lines to inter-cluster filaments traced by cluster-pairs. 

In this paper we use the redshift number density of absorp¬ 
tion lines, dN/dz, as a function of cluster-pair separation, as the 
relevant statistical quantity to characterize inter-cluster filaments 
(if any). We have chosen dN/dz as opposed to the number of 
systems per absorption distance, dN/dX, (or both), only for 
simplicity. Still, in Appendixj^we provide tables with relevant 
quantities and results for both dN/dz and dN/dX. We note that 
our conclusions are independent of this choice. 

At this point, it is also important to emphasize that we do 
not know a priori that cluster-pairs in our sample are tracing true 
inter-cluster filaments, and that even if they do, we do not know 
if these could produce a signal in the observed incidences of H I 
and O VI absorption lines at the S/N level obtained in our Q1410 
HST/COS spectrum. Although cosmological hydrodynamical 
simulations suggest that this may be the case, this paper aims to 
provide a direct test of such an hypothesis. Therefore, we will 
explore the behaviour of dN/dz over a wide range of scales 
around cluster-pairs both along and transverse to the Q1410 LOS. 
This means that in the following, we will allow the maximum 
impact parameter for clusters and cluster-pair inter-cluster axes to 
the Q1410 sightline to be larger than the fiducial values adopted 
in Section [33 and Section [ 3 ^ (i.e. larger than 20 and 3 Mpc, 
respectively). 

5.1 Measuring the redshift number density of absorption 
lines around cluster-pairs 

We measure the redshift number density of absorption lines 
around cluster-pairs in the following way. Let Ad be the maxi¬ 
mum impact parameter between a cluster-pair inter-cluster axis 
and the Q1410 sightline in a given calculation. Let Av be the 
maximum rest-frame velocity window to a given cluster-pair, at 
the redshift of such cluster-pair (see Figure[T]for an schematic 
diagram). Then, we define N{Ad, Av, lU™'") as the number of 
absorption components found within Av from the closest (in rest- 
frame velocity space) cluster-pair, from those cluster-pairs being 
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Full sightline 



Figure 5. Redshift path, Az, as a function of minimum rest-frame equivalent width for our survey of H I (black solid line and circles) and O VI 

(green dashed line and squares) absorption lines. The left panel shows the corresponding redshift path associated to regions of our Q1410 HST/COS 
spectrum within rest-frame velocity differences An = 1000 kms“^ from cluster-pairs at impact parameters smaller than Ad = 3 Mpc. The right panel 
shows the total corresponding redshift path for the full Q1410 HST/COS spectrum between 0.1 ^ z Si 0.5. Vertical lines show the minimum rest-frame 
equivalent width for our detected H I (solid black; Wr = 0.039A) and O VI (dashed green; Wr = O.OOlA) absorption line samples, excluding those in 
the category ‘c’ (uncertain; see Section[43). See Section[5T]for further details. 


within Ad from the Q1410 sightline, having rest-frame equiva¬ 
lent widths Wr ^ VF™‘"P^Let Az{Ad, Av, be the red¬ 

shift path in which a given absorption line having Wr ^ 
could have been detected along portions of the spectrum being 
within Av to any cluster-pair, from those cluster-pairs being 
within Ad from the Q1410 sightline, and being within our red¬ 
shift range constraint (i.e. 0.1 ^ z ^ 0.5)P^Then, the redshift 
number density is calculated as. 


dN 

dz 


{Ad, Av, wr") = 


N{Ad, Av, Wr" 


Az(Ad, At), ■ 


(3) 


Our methodology for estimating Az{Ad, Av, IF™'") is pre¬ 
sented in Appendix [C] 

Figure shows the redshift path, Az, along the Q1410 
sightline as a function of minimum rest-frame equivalent width 
,1F™‘'', for our survey of H I (black solid line and circles) and 
O VI (green dashed line and squares) absorption lines. The left 
panel shows the corresponding redshift path associated with 
regions of our Q1410 HST/COS spectrum within rest-frame 
velocity differences Aw = 1000 km s“^ from cluster-pairs at 
impact parameters smaller than Ad = 3 Mpc, while the right 
panel shows the total corresponding redshift path for the full 
Q1410 HST/COS spectrum between 0.1 ^ z 0.5. We see that 
the completeness level is very similar between the portions of 
the spectrum close to cluster-pairs and that of the full spectrum. 
We checked that this is also the case for multiple choices of Aw 
and Ad values, increasing from our fiducial values to cover the 
full spectrum (not shown). The vertical lines in Figure|^show 
the minimum rest-frame equivalent width for our H I (black 
solid; wr" = 0.039A) and O VI (dashed green; 1F”‘" = 
O.OOlA) absorption line samples, excluding those labelled as 
‘uncertain’ (category ‘c’; see Section [43| l. We see that these 
values correspond to high completeness levels and are therefore 


In our analysis we will also impose the line to be detected at 
WrjSWr ^ 3 (see Section]^, but this is not a requirement of the 
methodology. 

This redshift range condition is not a requirement of the methodology 
in the most general case. 


adopted as the minimum equivalent widths in the forthcoming 
analysis (but note that we could have detected O VI down to 
Wr ~ O.O 3 A with a similar completeness). 

5.2 Estimating the field redshift number density from the 
Q1410 sightline 

In this paper we have introduced slightly different ways to count 
and assess the statistical significance of absorption lines com¬ 
pared with has been done in previous published works. This is 
so because we opted to do a uniform analysis for H I (either 
total, broad or narrow) and O VI absorption lines, while previous 
works have usually focused on one type at a time. Therefore, 
we estimate the held redshift number density of a given species 
using our own methodology using the Q1410 sightline data alone. 
This is justihed by the fact that cluster-pair hlaments (if any) 
should only influence specific portions of the spectrum (in our 
case about ~ 1/6 of it), while the rest should match the held 
expectation (i.e. that from a randomly selected sightline). 

As described in section Section [Td] our sightline is ex¬ 
tremely unusual in terms of the number of cluster-pairs close to it 
(by construction, see also Section [2.2.1[ l. Therefore, an estimation 
of the held value from this sightline alone could be biased when 
compared against a representative ensemble of random sightlines. 
In order to correct for this potential bias we have proceeded 
as follows. Let Ntot be the total number of relevant absorption 
lines in the full Q1410 sightline between 0.1 ^ z ^ 0.5, and 
Wcpairs be the number of such absorption lines associated with 
our cluster-pairs assuming fiducial values of Aw = 1000 km s“^ 
and Ad = 3 Mpc (see Section|5T]for a definition of these two 
quantities). Therefore, our expected number of absorption lines 
associated with the held value can be estimated by, 

( Mcpitr \ 

Wfleld ~ Wtot — Wcpairs ^1-<51410 1 ’ 

where nj^a1r° number of independent cluster- 

pairs found in Q1410 sightline and those randomly expected, 
respectively. In our case = 7 and nfpau ~ 1 ± 1. 

Therefore, we take the expected field value to be Afield ~ 
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Impact parameter. Ad (Mpc) 




Figure 6. Redshift number density of total H I (brown circles; top panels) and 0 VI (green squares; bottom panels) absorption components as a function of 
rest-frame velocity window (Av; left panels) and maximum impact parameter to the closest cluster-pair axis {Ad', right panels) for a fixed Ad = 3 Mpc 
and Av = 1000km s~^, respectively. Note that bins are not independent from each other, as emphasized by the coloured areas. The total number of 
lines and redshift paths per bin are given in grey numbers on top of the datapoints. The expected field value estimated from our Q1410 sightline is 
represented by the horizontal dashed line with its ±lcr uncertainty represented by the darker grey region. The darkest grey hashed regions represents the 
±lcr field values from the|Danforth & ShullH2008t survey. See Section[6T]for further details. 


Ntot - (0.85 ± 0.14) Vcpairs[3 Finally, we estimate the rel- Q1410 alone matches those of comparable previously published 

evant redshift number density by using this corrected Vfleid as, blind surveys, making our assumptions reasonable. 


dN Vfleld 


(5) 


where Az is the total redshift path associated with the full Q1410 
sightline between 0.1 ^ 2 Si 0.5. 

This methodology assumes (i) that there is an excess of ab¬ 
sorption lines in the data compared to the field expectation, and 
(ii) that this excess is purely confined within An = 1000 km 
and Ad = 3 Mpc from the known cluster-pairs. If assumption 
(i) is incorrect, then our field expectation estimation will be un¬ 
derestimated. If assumption (i) is correct, but assumption (ii) is 
incorrect, then our field expectation estimation will be overesti¬ 
mated. In Sectionl^we show that our field estimations based on 


5.3 Statistical uncertainty estimations 


The statistical uncertainty in our calculations is dominated by 
the uncertainty in N{Ad, Av), which we assume is Poissonian 
and est imate from the analytical approximation g iven by|Gehrels| 
■jf « y/N + 3/4 -I- 1 and cr)^ « y/N — 1/4. The 


1 1986b: (7+ 


statistical uncertainty in our estimation of Afleid is taken from 
the contributions of both the Poissonian uncertainty of Ntot, 
and the statistical uncertainty of which we propagate 

assuming independence between these two quantities. Given that 
the statistical uncertainties in Az{Ad, Av) and AX (Ad, Av) 
are much smaller, we neglect them. 


We note that a consistent correction factor of Nftoid = Ntot — 
(0.85 ± 0.2) A^cpairs found when considering the total number 
of cluster-pairs instead of independent ones, i.e. = 11 and 


6 RESULTS 

In this section we report our results on dN/dz, for our different 
samples of H I (total, narrow and broad) and O VI absorption 
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Figure 7. Same as Figure|^but for narrow (6 < 50 km s“l; blue triangles; top panels) and broad (6 ^ 50 km s“l; red pentagons; bottom panels). The 
darkest grey hashed regions in the bottom panels represents the itlcr field value from the |Danforth et al.|(2010) BLA survey. See Section [6T] for further 
details. 


lines observed in the Q1410 sightline (see Section]^ applying 
the methodology described in Sectionj^to associate them with 
cluster-pairs (see Section]^. For simplicity, and in order to reduce 
the ‘shot noise’ of the measurements, the following results are 
obtained by varying An and Ad for fixed values of Ad and An, 
respectively (as opposed to varying both values at the same time). 

A summary of all the results presented in this section (and 
those of dN/dX, not described here), are given in Tables [ gT] to 
lG4l 


6.1 Redshift number densities 

Figure]^ shows the dN/dz of total H I (top panels; brown cir¬ 
cles) and O VI (bottom panels; green squares) absorption com¬ 
ponents as a function of maximum velocity window (An; left 
panels) and maximum impact parameter (Ad; right panels) for 
a fixed Ad = 3 Mpc and An = 1000 km s“^, respectively. The 
expected field values following our approach described in Sec¬ 
tion |5^ are shown by the horizontal dashed line with its Tier 
uncertainty represented by the grey region. We also show the 
±1(T field values from the |Danforth & ShuIT| ( |2008| l survey as the 
darker grey hashed regions, which is consistent with ours. 

When we fix Ad = 3 Mpc (left panels), we observe a 


clear overall increase in the redshift number density of H 1 and 
O VI absorption lines with decreasing An. Similarly, when we 
fix An = 1000 km s“^ (right panels), we observe an overall 
increase in the redshift number density of H 1 and O VI absorption 
lines with decreasing Ad, but only down to Ad ~ 3 Mpc; at 
Ad < 3 Mpc a flattening (or even decrease) trend is observed, 
which we believe is mostly due to our small sample in such 
bins|^This change of behaviour motivated our adopted fiducial 
value of 3 Mpc for the maximum transverse separation between 
cluster-pair axis and the Q1410 sightline (see Section [T2l >. 

To test whether kinematic trends are present in the H 1 data, 
we repeated the dN/dz measurements for both narrow (NLA; 
h < 50km s“^) and broad (BLA; b 50km s“^) H I Lya 
absorption lines. Although the canonical value for BLAs tracing 
the WHIM is 40 km s“^, this limit assumes that the broadening 
is purely thermal. Following more recent work ( [Richter et al.| 
|2006[|Danforth et al.|2010| l, it is acknowledged that non-thermal 
broadening mechanisms are likely to be present in absorption 
line samples. Thus, our adopted value of 6 ^ 50 km s“^ is more 
conservative (see also Section[7^. 

But note that with this limited sample we cannot rule out that a real 
decreasing signal is present either. 
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Figure 8. Relative excesses of redshift number densities compared to the field expectation, as a function of rest-frame velocity window (An; left panels) 
and maximum impact parameter to the closest cluster-pair axis (Ad; right panels), for fixed Ad = 3 Mpc and An = 1000 km s“ ^ values respectively. 
The top panels show those for total H I (brown circles, solid lines) and O VI (green squares, dashed lines) absorption line samples. The bottom panels 
show those for narrow (NLA; b < 50 km s“^; blue triangles, solid lines) and broad (BLAs; b ^ 50 km s“^; red pentagons, dashed lines) H I absorption 
line samples. Coloured light shaded areas represent ztlc statistical uncertainties. See Section[6!^for further details. 


Figure|^is equivalent to Figure|^but for NLA (top panels) 
and BLA (bottom panels) absorption line samples. The hashed 
darker grey area in the bottom panels represents the field value 
obtained by |Danforth et al.|(2010| > for BLAs. 

When we fix Ad = 3 Mpc (left panels), we observe a clear 
overall increase in the redshift number density of both narrow 
and broad FI I absorption lines with decreasing An. When we fix 
An = 1000km s“^ (right panels), we also observe an overall 
increase down to Ad ~ 3 Mpc; below this scale a decreasing 
trend may be present for narrow FI I lines, while for broad FI I 
lines the increasing trend persists. 


We also note that our estimation field expectations are fully 


consistent with those from previous 

blind surveys (|Danforth| 

& Shull||2008 

Danforth et al. 

2010 l 


This implies that our 


characterization of absorption lines (see Section and our 
methodology for estimating the field expectation from our 
Q1410 data alone (see Section [5^ are reasonable. Therefore, 


Note that in the case of BLAs, both field values have similar un¬ 
certainty. This is because PPanforth et al.|j2()10^ included a systematic 
contribution to the error, whereas ours is purely statistical. 


we can conclude that the vast majority (if not all) of the observed 
excesses come from scales within An = 1000 km s“^ and 
Ad = 3 Mpc (see Section|7]for further discussion). 


6.2 Relative excesses with respect to the field 

Figure[^shows the relative excesses of redshift number densities 
of our absorption line samples compared to their respective field 
expectations, defined as ^ ^ | field, as a function of rest-frame 
velocity window (An; left panels) and maximum impact param¬ 
eter to the closest cluster-pair axis (Ad; right panels), for fixed 
Ad = 3 Mpc, An = 1000 km s“^, respectively. The top panels 
show the results for our total H I (brown circles, solid line) and 
O VI (green squares, dashed line) samples, while the bottom 
panels show the results for our NLA (blue triangles, solid line) 
and BLA (red pentagons, dashed line) samples. Coloured light 
shaded areas represent the Tier statistical uncertainties. 

Although subject to large statistical uncertainties, the rel¬ 
ative excess for BLAs tends to be the highest of all, reaching 
a value of ~ 6 times its field expectation at Ad = 3 Mpc and 
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Figure 9. Representation of the position of the Q1410 sightline with 
respect to different independent cluster-pairs in our sample. The y-axes 
conespond to the impact parameter to a given cluster-pair, Ad, while the 
x-axes correspond to the distance to the closest galaxy cluster in a given 
cluster-pair, along the cluster-pair axis, Aa;, (see Figure[T]for a schematic 
diagram). The four panels correspond to different absorption line samples, 
from top to bottom: total H I narrow H I (NLAs; b < 50 km s”^), 
broad H I (BLAs; b ^ 50 km s“^), and total O VI. Filled symbols 
correspond to portions of the sightline showing absorption lines within 
An = 1000 km s“^ from the redshifts of a given cluster-pair, where 
the sizes of the symbols are proportional to the number of absorption 
lines. Empty symbols correspond to portions of the sightline showing 
no absorption. When multiple cluster-pairs lie at similar redshifts (i.e. 
from grouped ones), the absorption lines are associated to the one having 
the smallest impact parameter value and the rest of the cluster-pairs are 
obviated. We also show the typical virial radii of our sample of clusters 
with a quarter dashed circumference of radii 1 Mpc centred at the origin. 
See Section lTTI for further details. 


Av = 1000 km s“^. On the other hand, the excess of NLAs 
tends to be the smallest of all, reaching a value of only ~ 1.5 
times its field expectation at the same scales. The relative excess 
of total H I tends to lie in between that of NLAs and BLAs, but is 
closer to that of NLAs because these type of absorbers dominate 
the neutral hydrogen sample. The sample of O VI has the largest 
statistical uncertainties (it is indeed the smallest sample), which 
makes its relative excess to be consistent with all others even at 
the 1(7 c.l. Strictly, within ~ 2a c.l. all the reported excesses are 
consistent with each other across different samples, and are also 
consistent with their respective field expectations. Therefore, it 
is important to test these trends with larger datasets. 


6.3 Equivalent widths distributions 

In Appendix [D| we provide a comparison between equivalent 
widths distributions for our different samples. We did not find sta¬ 
tistically significant differences between systems close to cluster- 


pairs and the field expectation in terms of equivalent widths, at 
least from our limited sample sizes. 


7 DISCUSSION 

7.1 Filamentary structure 

Here we argue that our results are roughly consistent with a fila¬ 
mentary structure for the absorbing gas close to cluster-pairs. 
This is so because when we restrict the analysis to a fixed 
Av = 1000km s“^, the excess is maximized at impact parame¬ 
ters of Ad ~ 3 Mpc (or even Ad < 3 Mpc for broad H I), while 
the typical separation between clusters in our cluster-pairs are of 
the order of > 10 — 15 Mpc (see sixth and seventh column of 
Table|^. Moreover, we have also found that when we consider 
scales far outside our fiducial values, we fully recover the field 
expectation (see reasoning presented in Section|5^. From the cu¬ 
mulative results presented in Section]^ we can directly calculate 
dN/dz in independent intervals instead by subtracting both the 
reported number of absorption lines and redshift path in a given 
bin (i.e. Wbin and Azbin) of a smaller scale, to those at the scale 
of interest. For instance, if we focus on the total H I sample and 
consider scales between 3 to 100 Mpc as those of interest, we can 
estimate dN/dz as ^ ~ (o^257o^o6) ~ ~ 116 ± 25 (see 

top panel of Figurej^and Table [GT| . This number is fully con¬ 
sistent with the field expectation (e.g. |Danforth & Shull|200^ , 
and therefore we conclude the vast majority (if not all) of the 
observed excesses come from scales within Av = 1000 km s“^ 
and Ad = 3 Mpc. However, because we did not impose a mini¬ 
mum distance between the QSO sightline and the closest cluster 
of a cluster-pair, there is also the possibility that our survey ge¬ 
ometry does not represent that of a filamentary structure. If those 
separations are all < 3 Mpc for instance, our survey could be 
probing a more spherical (or disk) geometry instead. 

Figurej^shows a geometrical representation of our survey. 
We observe that the distances probed by our survey cover scales 
between Aa; ~ 0 — 10 Mpc along the inter-cluster axes, roughly 
uniformly. Therefore, we conclude that the geometry of our 
survey is indeed consistent with that of a filamentary structure, 
but we also stress that a larger sample must be analysed in order 
to better constrain the geometry. 


7.2 Covering fractions of absorbing gas close to 
cluster-pairs 

Here we provide a first estimation of the covering fractions, 
of the absorbing gas close to cluster-pairs and compare 
them with the random expectation, Our adopted fiducial 

Av = ±1000 km s“^ corresponds to ~ ±16 Mpc along the 
line-of-sight (if cosmological). This is a larger scale compared 
to our fiducial filament radius of ~ 3 Mpc. Therefore, the ex¬ 
cesses do not necessarily come from single inter-cluster filaments. 
Although our reported dN/dz signals tend to keep increasing 
at smaller rest-frame velocity differences, the samples also get 
smaller, which makes the statistical uncertainties larger too (see 
Figures|^and|7]l. By comparing the observed covering fractions 
to random expectations, we can shed light into the origin of the 
reported excesses in relation to the cluster-pairs themselves. 

From Figure]^ we observe that 7/7 sightlines close to 
cluster-pairs did show at least 1 H I absorber, which implies a 
covering fraction of /f“‘(HI) ~ 1.00tg;^!^|^Similarly, NLAs 
and BLAs were both found in 5/7 of the sighthnes probing them 


The uncertainty is estimated assuming a binomial distribution for 
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Figure 10. Estimated covering fractions of gas close to cluster-pairs 
solid symbols) and in the random expectation open 

symbols), for our different samples of absorption lines using the same 
symbol/colour convention as in previous figures. We observe that broad 
H I absorbers (6 ^ 50 km s“^) have about ~ 4 times a larger covering 
fraction in close to cluster-pairs than from the random expectation. See 
Section[T2]for further details. 

(althou^ different subsamples; see second and third panel of 
Figureg, implying /f“‘(NLA) « /«”^(BLA) r. 0.7llg;«. 
In contrast, O VI absorbers were found in 1/7 of the sight¬ 
lines probing them, implying a smaller covering fraction of 
jflmt(ovi) ~ 0.141 q ^ 2 - These results are summarized in Ta- 
ble|^(upper half) and shown in Figure p^(filled symbols), using 
the same symbol/colour convention as in previous figures. 

As a control sample for a given absorber, we estimated 
as. 


A2 ;(Av) ^Ifieid if Az(At;)^|fleid < 1 ,,, 

1 if Az(Aw)^|fleid > 1 

where ^ | fleid is the redshift number density of lines in the field, 
and A 2 (Ar;) is the corresponding redshift path to a rest-frame 
velocity window Av evaluated at 2 : = 0.35 (the median redshift 
of our cluster-pair sample). We note that we are neglecting the 
intrinsic clustering of absorption lines with themselves, which is 
justified because our estimations are obtained from inde¬ 
pendent structures]^ The fifth column of Table[^(upper half) 
summarizes the random expectations for our different samples 
using Av = dzlOOO = 2000km s“^. These results are also 
shown in Figure|T^as open symbols. We observe that the cover¬ 
ing fractions for total H I, NLAs and O VI close to cluster-pairs 

the number of ‘hits’ given the 7 independent ‘trials’, using the Bayesian 
formalism described by|Cameron|^201 1| with a flat prior. 

For obtained from non-independent structures (e.g. a single 

well mapped filament), a meaningful estimation of ^nust take 

clustering into account. At first order, there will be two opposing effects 
(when clustering is positive): clustering along the line-of-sight will tend 
to decrease while the clustering transverse to the line-of-sight 

will tend to increase it. Moreover, higher order correlations must also be 
considered to account for the joint probability of having multiple ‘hits’ 
in such large single non-independent volume. 


Table 5. Estimated covering fractions of absorbing gas close to 
cluster-pairs. 


Sample 

(1) 

^hits 

(2) 

"^trials 

(3) 

y?flmnt 

(4) 

yrand 

(5) 

Excess 

(6) 


Av = 

±1000 km s ^ and Ad 

= 3Mpc 



Hi 

1 

7 

1 OO"*”^'®^ 

i.UU_Q 23 

1-UU_Q Q5 


1 

NLA 

5 

7 

Q yi + O.lS 
' ^- 0.26 

0 qo + 0.07 


1 

BLA 

5 

7 

0 71 
^ -*-- 0.26 

0 1 fi +010 

U.iD_o 10 


4 

O VI 

1 

7 

^•^^-0.12 

0 1 5 + 0.03 
^•-*-^- 0.03 


1 


Av = 

: ±500 km ; 

and Ad = 

= 3 Mpc 



Hi 

5 

7 

0 71 + 0-18 

O'-' -*-- 0.26 

0 58+0 03 
^•^°- 0.03 


1 

NLA 

3 

7 

n 4Q+0-25 

o o 
6 d 
+ 1 

o 


1 

BLA 

4 

7 

0 57 + 0-22 
- 0.25 

0.08lg;°5 


7 

0 VI 

1 

7 

0 1 4 + 0-26 

0-081^2 


1 


(1) Sample of absorbing gas. (2) Number of ‘hits’ defined as 
sightlines showing absorption in a given sample, within Av = 
{dilOOO, ±500} km and within Ad = 3Mpc. (3) Number 
of ‘trials’ defined as the total number of sightlines to look for 
absorption. (4) Covering fraction close to cluster-pairs estimated 
as T^hits/^trials (uncertainties correspond to those of a binomial 
Icr C.I.). (5) Covering fraction in a random sightline for a given 
Av = {±1000, ±500} km s~^. (6) Excess covering fraction de- 
fined as 


are consistent with their random expectations, while //*™'^(BLA) 
is about ~ 4 times larger than its random expectation (at the 
~ 2a C.I.). Having consistency with the random expectations are 
not surprising; as mentioned, Av = ±1000km s“^ corresponds 
to about ~ ±16Mpc along the line-of-sight (if cosmological) 
around massive structures traced by galaxy cluster pairs. In this 
scenario, we do expect that the dN/dz\fieid values to be domi¬ 
nated by absorption lines found in the overdense LSS. On the 
other hand, having an excess in the covering fractions of BLAs 
with respect to the random expectation implies that this type of 
gas is not common over ~ ±16Mpc scales around cluster-pairs, 
and therefore it has to come from smaller scales (i.e the cluster- 
pair itself). To test this hypothesis, we have repeated the covering 
fractions estimations using a smaller An = ±500 km s“^, and 
the excess of BLAs remains large (~ 7; see values in the bottom 
halfof Tablel^. This behaviour favours the conclusion that most 
of the BLAs in our sample are directly related to the cluster-pairs 
themselves. 


7.3 Could the observed excesses of gas be due to galaxy 
clusters/groups or individual galaxy halos? 

Regarding massive structures, we can see from Figure|^that the 
vast majority of sightlines are probing regions far away from the 
virial radii of known galaxy clusters (see also the left panels of 
Figure [FT| . However, because of the limited completeness level 
of the redMaPPer catalog, there could still be unknown clusters or 
groups in such regions. In order to directly address this question, 
one must survey the Q1410 field for individual galaxies and LSS 
close to the Q1410 line-of-sight, which we leave for future work 
(see Section|T^. Still, in the following we provide an assessment 
of the plausible incidence of gas associated to individual halos 
from two indirect but independent arguments. 
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7.3.1 Redshift number density of galaxies 


Based on the reasoning presented by |Prochaska et ar] ( |201 l| l, 
one can estimate the redshift number density of galaxies of 
luminosity L Js Lmin, dN/ dz\g^z.iB, intersecting a given LOS 
by assuming a given cross-section for them. For galaxies with 
L Js 0.01 — O.OOIL* and assuming cross-sections given ^ 
the virial radius of galaxies with unity covering fractions!^ 
Prochaska et al. i2011 i find that dN/dz\^BiB ~ 10 — 20 (see 


their figure 8). This estimate is valid for the field and so we 
need to correct for the fact that cluster-pairs generally probe LSS 
overdensities. We use the excess of narrow H I lines to estimate 
the overall overdensity in our cluster-pair sample, as ~ 1 — 2 
times the mean density of the Universe (note that these values 
are consistent with the expected overdensities traced by cluster- 
pairs from cosmological simulations). Therefore, we estimate 
dN/dz\lX^ ~ 10 — 40. This number is lower than the typical 
dN/dz ~ 60 — 300 observed for our samples of H I (either 
narrow or broad), but is comparable to the dN/dz ~ 30 — 100 
observed for O VI. Therefore, it seems unlikely that most of the 
excess observed for narrow and broad H I gas close to cluster- 
pairs is driven by galaxy halos of individual galaxies. On the other 
hand, our reported excess of O VI gas might well be produced 
(at least partly) by individual galaxy halos (see also below). 


7.3.2 Metal absorption lines 

One can also infer the presence of galaxy halo material by 
means of metal absorption lines, in particular from low-ionization 
species (e.g. |Werk et al.|20T4l l. There is only 1 absorption system 
in our cluster-pair sample showing metal absorption lines: the 
one at z « 0.35 from which all the 3 observed O VI components 
come from (see the fourth panel in Figure [FT| . This system has 
strong FI I absorption with column densities N > 10^'‘cm"2 
and shows a complex kinematic structure (4 components in total, 
3 narrow and 1 broad). The second narrow component also shows 
the presence of C II, Si III and N III (see Table]^. Therefore, 
although it seems very likely that an important fraction of the 
absorbing gas in this system comes from an individual galaxy 
halo or its immediate surroundings, this system only accounts for 
3/10 ~ 30 per cent of NLAs and 1/6 ~ 16 per cent of BLAs in 
our sample. Again, from this independent reasoning we reached 
the same conclusion as before, i.e. individual galaxy halos could 
account for the observed excess in O VI lines (although not a 
requirement), but not for the majority of H I gas. 


7.4 Statistical evidence of the WHIM 

Here we consider whether the observed trends could also be 
consistent with the presence of a warm-hot intergalactic medium 
(WHIM). Ideally, one would require a full characterization of 
the physical conditions of individual absorbers using multiple 
species and comparing them with the expectations of different 
models of ionization. However, because our sample is dominated 
by absorption systems having no other species than H I, this 
approach is not feasible. Even when other species are present in 
individual systems (e.g. H I and O VI) this approach requires the 
uncertain assumption that the majority of the gas comes from a 
single phase, which is controversial at the very least (e.g. Werk et 
al., in prep.). In view of this intrinsic limitation for an individual 
characterization, here we opted for a purely statistical approach. 

The WHIM is usually defined as gas at temperatures in the 


range of T ~ 10® — 10^ K, implying a minimum Doppler param¬ 
eter of 6 = 40km s“^ for H I (i.e. assuming the broadening is 
purely thermal). Therefore, H I lines with h < 40 km s“^ can not 
be caused by a WHIM. Non-thermal processes can also broaden 
absorption line profiles, including turbulence, Hubble broadening 
and unresolved blends (e.g. |Garzilli et al.|2015[ and references 
therein). In overdense environments, we expect turbulence to be 
the dominant source of non-thermal physical broadening (but 
see Section [73] for a discussion regarding unresolved blends). 
Assuming that the turbulence contribution is higher for hotter 
gas, we should have feturb ~ a&T, and hydrodynamical sim¬ 
ulations suggest 0 ^ a ^ 1 (e.g. [Tepper-Garcfa et al.|2012| l. 
Assuming the most extreme case of a = 1, we have that the ob¬ 
served Doppler parameter would be & = s/2bT, making a limit 
b > 40\/2 « 57km s“^ extremely conservative in ensuring to 
trace HI gas at T > 10® K. However, there could still be genuine 
WHIM H I absorption lines in the 40 < fenm < 40\/2km s“^ 
range (i.e. when thermal broadening does dominate). 

In this paper, we used a limiting Doppler parameter value 
of feiim = 50 km s“^ instead, which is in between 40 < feum < 
40\/2km s“^. This limit was partly chosen for allowing a di¬ 
rect comparison to previous published work (e.g.jDanforth et al.j 
|2010| ), but also because it minimizes potential misidentification 
of lines that are supposed to trace warm-hot gas but trace cool 
gas instead, and viceversa (i.e. genuine warm-hot absorbers hav¬ 
ing b < 6iim, and genuine cool absorbers having b ^ &lim)[f] 
Indeed, using foiim = 50km s“^, there is only 1 (1) line in the 
range 40 < b < 50km s”'^ (50 < b < 40x/2km s“^) in our 
NLA (BLA) sample associated with cluster-pairs (e.g. see points 
with white stars in Figure]^. A limit of biim = 50km s“^ also 
corresponds to a a = 0.75, and therefore is still quite conserva¬ 
tive in ensuring that BLAs trace gas at T 10® K, even with 
a substantial turbulence contribution. Regardless of these con¬ 
siderations, we also note that about half of the BLAs associated 
with cluster-pairs in our sample are actually very broad, with 
Doppler parameters 6 ~ 80 — 150km s“^ (e.g. see Figure]^, 
which should make them more likely to trace gas at WHIM 
temperatures (but see Section [73| >. 

One of our proposed diagnostics is to compare the excesses 
in the incidence of narrow and broad H I absorbers (and even¬ 
tually O VI when larger samples are gathered) found close to 
cluster-pairs with respect to their field expectations (see Figure]^. 
Because inter- cluster filaments correspond to overdense regions 
in the Universe, an excess of gas is generally expected to occur, 
and indeed we have shown that this is the case (see Section]^. 
Under the null-hypothesis that BLAs and NLAs probe gas in 
similar physical conditions (i.e. similar physical entities), then 
we expect both these excesses to behave in a consistent manner. 
On the other hand, if BLAs and NLAs are not probing similar 
physical conditions, a different behaviour for the excesses is ex¬ 
pected instead. A WHIM signature associated with inter-cluster 
filaments may include the relative excess of BLAs to be higher 
than that of NLAs, which is exactly what we observed (although 
only at the ~ Icr c.F; see Figure]^ bottom panels). By increas¬ 
ing the sample sizes we may test for any statistically significant 
difference between them. 

Another proposed diagnostic is to constrain the overall ge¬ 
ometry for the excess of gas, in terms of both BLAs and NLAs 
(and O VI when larger samples are gathered). A WHIM signature 
in this context would require an increase in the covering fractions 
of BLAs towards the inter-cluster axes compared to the random 
expectations, which is also what we have observed (again, only 


Which is a conservative assumption (e.g. seejWakker & Savage|2009| 
[Prochaska et al.|2011[|Johnson et al.|2014[ for counter examples). 


This is not the same as choosing the limit that maximizes the differ¬ 
ence between observed incidences though (which we did not try). 
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at the ~ 2(t c.L; see Section|T^. Moreover, one can also look 
for trends in the Doppler parameters of H I absorption lines with 
respect to impact parameter to the inter-cluster axes, as a proxy 
of temperature. Assuming simple models for the ionization of 
the gas (e.g. purely collisional), one can even use the inferred 
temperatures to estimate a total hydrogen column density from 
the observed H I one. A WHIM signature should produce, on 
average, higher hydrogen column densities for higher tempera¬ 
tures. 

Although all our tentative results (~ 1 — 20- c.l.) may be con¬ 
sistent with the presence of a WHIM in inter-cluster filaments, we 
emphasize that a larger sample must be analysed before reaching 
a definite conclusion. 

7.5 Caveat 

Probably the most important source of concern in our experi¬ 
mental design, is our limited ability to disentangle blends, which 
is key to detect broad and shallow absorption features expected 
to arise in the WHIM. The importance of this systematic uncer¬ 
tainty depends on the S/N, as the higher the S/N the easier it is 
to assess the kinematic structure of the absorption feature (e.g. 
[Richter et al.|2006{ [Danforth et al.|20I0| l. We emphasize that 
this is an intrinsic limitation of the absorption-line technique, 
meaning that all these kind of observational samples are, to some 
extent, affected by this issue. As described in Section|^ we at¬ 
tempted to avoid this bias by fitting asymmetric lines with at 
least two components. Although not impossible, we believe that 
the likelihood of having misidentified multiple narrower blended 
components as a single broader and symmetric one in a large 
fraction of our H I sample is low (see Figure [FT] for individual 
examples close to cluster-pairs). Given that we are comparing the 
relative incidences of lines between different samples drawn from 
spectra of similar S/N (see Section]^, our statistical approach 
seems adequate for minimizing this potential source of uncer¬ 
tainty (as opposed to attempting a full physical characterization 
of individual systems). 

7.6 Future prospects 

Despite our promising results, the existence of the WHIM in inter¬ 
cluster filaments still needs to be observationally confirmed; our 
pilot survey was not design to draw statistically significant results, 
but primarily to show that such a goal is currently possible with 
existing instrumentation. In this section, we enumerate remaining 
work for providing a firm detection of the elusive WHIM in the 
context of our methodology. 

7.6.1 Increase the sample sizes 

Increasing the sample sizes is a key requirement. In this re¬ 
spect, we are actively working on two fronts: (i) pursuing new 
HST/COS observations of targeted QSO sightlines intersecting 
multiple cluster-pairs; and (ii) searching in the HST archive for 
already observed QSOs whose sightline intersects single or mul¬ 
tiple cluster-pairs. We estimate that the HST/COS archive will 
allow us to extend the present work to tens of sightlines, but 
approach (i) will still be necessary for efficient follow up obser¬ 
vations (e.g. galaxy surveys; see Section [7.6.2| l. 

7.6.2 Survey for galaxies around the QSO sightline 

As discussed in Section [73| we need to survey galaxies around 
our QSO sightline in order to directly rule out the potential asso¬ 
ciation of BLAs with the halos of individual galaxies. To this end. 


spectroscopic redshifts are needed (current galaxy photometric 
redshift uncertainties are too large for meaningful associations 
with absorption lines). We are currently pursuing multi-object 
spectroscopy (MOS) and integral field units (IFU) observations 
around the Q1410 field. We will use MOS surveys to assess the 
distribution of galaxies over 0.3 — 10 Mpc scales. This will be 
important not only for determining whether galaxy groups or 
clusters are responsible of our observed excesses of absorbing 
gas, but also to determine the actual geometry of the LSSs inter¬ 
sected by the Q1410 sightline, including: (i) assess whether these 
putative inter-cluster filaments are straight or bent; and (ii) refine 
the cosmological redshift of the structures at the position of the 
Q1410 sightline. The IFU observations will primarily focus on 
mapping galaxies on scales within < 100 — 300 kpc to the QSO 
sightline (the typical CGM scales) at a very high completeness 
level, including faint star-forming galaxies with no detectable 
continuum but having bright enough emission lines. 


7.6.3 Comparison to hydrodynamical simulations 

Another key aspect of this project, is the comparison of our obser¬ 
vational results to those obtained from cosmological simulations. 
Our experimental design offers a unique opportunity to test the 
prediction of the ACDM paradigm in the largest and densest 
filaments of the cosmic web, while also constraining current 
galaxy evolution models. As we have shown, H I dominates the 
gaseous content found in inter-cluster filaments (see Section]^, 
and is very likely that they originate far away from individual 
galaxy halos (see Section |7.3| ). This makes a direct compari¬ 
son to simulations straightforward because we expect this gas 
to be unaffected by the uncertain baryonic processes occurring 
in galaxies (e.g. SNe/AGN feedback). On the other hand, the 
subsample of absorption systems showing metal absorption (e.g. 
those with O VI) can put constraints on these uncertain ‘sub-grid 
physic’ models for galaxy formation. However, the need for a full 
treatment of shocks in the gas limits the approach to being hydro- 
dynamical. If the predictions of WHIM signatures in inter-cluster 
filaments match our observational results, they will provide yet 
another piece of evidence supporting the existence of this elusive 
medium. 


8 SUMMARY 

The warm-hot intergalactic medium (WHIM) has been predicted 
to account for a significant fraction of the baryons at low-z, but 
its actual existence has eluded a firm observational confirmation. 
In this paper, we have presented a novel approach for detecting 
the WHIM, by targeting regions of the cosmic web where its 
presence is predicted to be ubiquitous: the putative filaments 
connecting galaxy clusters. As a proof of concept, we selected a 
single bright QSO (namely Q1410), whose exceptional sightline 
passes within Ad = 3 Mpc from the 7 inter-cluster axes con¬ 
necting 7 independent cluster-pairs at redshifts 0.1 ^ z 0.5, 
and observed it with HST/COS. We performed a full character¬ 
ization of absorption features in the FUV spectrum of Q1410, 
independently of the presence of known LSS traced by the galaxy 
cluster-pairs. From this dataset, we conducted a survey of dif¬ 
fuse gas along the QSO sightline with special focus on H I and 
O VI absorption lines. This survey allowed us to provide, for 
the first time, a systematic and statistical measurement of the 
incidence, dN/dz, of intervening H I and O VI absorption lines 
close to cluster-pairs. We split the sample of H I Lya into broad 
(BLA) and narrow (NLA) using a Doppler parameter limit of 
feiim = 50 km s“^, which ensures BLAs to be mostly from gas 
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at temperatures T ^ 10® K, even when accounting for turbu¬ 
lence. We quantified the incidence of H I, NLAs, BLAs and 
O VI absorption lines close to cluster-pairs by varying the mini¬ 
mum rest-frame velocity window, Aw, and the minimum impact 
parameter to the inter-cluster axes, Ad, and found that the inci¬ 
dence of diffuse gas is maximized at An < 1000 km s“^ and 
Ad < 3 Mpc. At these scales we report: 

• dA/d 0 (HI)|flnint = 287'liYi, which corresponds to 

2 .slo g times the field expectation; 

• dAt/d 2 (NLA)|flinnt = ITOlgg which corresponds to 
l-7lo 6 times the field expectation; 

• dAt/d 2 (BLA)|flmnt = IO 8 I 43 which corresponds to 
B.llg'j times the field expectation; and, 

. dA/d 0 (OVI)|fl„nt = 581” which corresponds to 

4 .ll 2 y times the field expectation. 

Although individually these excesses are only at the ~ 
1 — 2a C.I., in concert they suggest a physical overdensity 
close to cluster-pairs. Our results are also roughly consistent 
with a filamentary geometry for the absorbing gas connecting 
cluster-pairs, with covering fractions of:/c(HI)Iflmnt = I.OIq'^; 
/c(NLA)|fl^„t = 0.7lgj; /c(BLA)|fl„„t = 0.7l°j; and 
/c(OVI)|fimnt = O.ldlg ®. Our resported covering fractions 
of total H I, NLAs and O VI are all consistent with their ran¬ 
dom expectations. In contrast, the reported covering fraction of 
BLAs is a factor of ~ 4 larger than the random expectation. 
Because a rest-frame velocity window of Av « ±1000 km s“^ 
corresponds to a rather large co-moving distance along the LOS 
(i.e ~ ±16 Mpc), having consistency with the random expecta¬ 
tion is not surprising, and suggests that most of the excesses of 
NLAs and O VI absorption come from the overall LSS overden¬ 
sities around massive structures traced by galaxy cluster pairs. 
On the other hand, a higher covering fraction of BLAs close to 
cluster-pairs compared to the random expectation, suggests that 
the excess of BLAs is physically associated to the cluster-pairs 
themselves. Indeed, we also reported covering fractions using a 
Av ~ ±500 km s“^ and reached the same conclusion. 

Based on statistical arguments, we also concluded that most 
of the reported excesses of NLAs and BLAs cannot be attributed 
to individual galaxy halos but rather to truly intergalactic material. 
In contrast, the reported excess of O VI most likely comes from 
gas associated to individual galaxy halos or their immediate 
surroundings. 

We argued that a behaviour in which BLAs show larger 
relative excesses compared to that of NLAs (as tentatively re¬ 
ported here), may be a direct signature of the WHIM, especially 
if identified in the regions close to cluster-pairs. With a larger 
sample of QSOs and a careful accounting of systematic effects, 
the technique we have presented here should therefore enable a 
firm detection of the WHIM in inter-cluster filaments. 
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APPENDIX A: PROPERTIES OF CLUSTERS CLOSE 
TO QI410 COMPARED TO A CONTROL SAMPLE 

In this section we investigate whether our subsample of clusters 
close to the Q1410 sightline is a fair representation of the cluster 
population in the whole redMaPPer catalog. To do so, we will 
compare the mass (estimated from Equation (|T]l) and redshift 
distributions of our sample to appropriate control samples drawn 
from the redMaPPer catalog. 

The left panel of Figure [AT| shows the mass distribution for 
our subsample of 57 clusters within 20Mpc of Q1410, having 
richness values > 10 and lying between 0.1 ^ 2 : 0.5 (thick 

black line), and the normalized mass distribution from a control 
sample satisfying only the richness and redshift range criteria 
(162 144 clusters in total; thin red line). We see no apparent dif¬ 
ference between these two distributions, and no statistically sig¬ 
nificant differences are detected either: the Kolmogorov-Smimov 
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Redshift 


Figure Al. Comparison between mass (left panel) and redshift (right panel) distributions between our subsample of 57 redMaPPer clusters within 20 Mpc 
from the Q1410 sightline, having richness values ^ 10 and lying between 0.1 ^ z 0.5 (thick black line), and normalized control samples satisfying 
only the richness and redshift range criteria (thin red line; see Appendix|^for details). A Kolmogorov-Smimov test to the unbinned distributions give no 
statistically signihcant differences between the samples. 


test over the full unbinned samples gives a ~ 60% probability of 
both being drawn from the same parent distribution. 

In the case of the redshift distribution, we can not directly 
compare that of the 162 144 clusters satisfying only the richness 
and redshift range criteria, to our subsample. This is so because 
our subsample of clusters are defined by a cylinder (i.e. constant 
volume as a function of redshift), rather than fixed solid angle 
(increasing volume as a function of redshift). In order to pro¬ 
vide an appropriate control sample, we looked at the redshift 
distribution of redMaPPer clusters having richness values > 10 
and lying between 0.1 ^ z 0.5, in 5000 cylinders of radius 
20 Mpc, selected at random positions between R.A. £ [140, 222] 
degrees and Dec. £ [4, 56] degrees (i.e . we ll within the SDSS 
footprint)]^ The right panel of Figure Al shows the redshift 
distribution for our subsample of 57 clusters (thick black line), 
and the normalized redshift distribution from our aforementioned 
control sample (thin red line). Again, no statistically significant 
difference is detected: the Kolmogorov-Smirnov test over the 
full unbinned samples gives a « 20% probability of both being 
drawn from the same parent distribution. 

From these comparisons we conclude that no strong bias is 
present in our subsample of clusters close to the Q1410 sightline. 


APPENDIX B: COMPARISON BETWEEN 
ABSORPTION LINE SIGNIFICANCE ESTIMATION 
METHODS 


Here we compare two methods for estimating the significance 
level of absorption features, namely our Wi / 5Wr criterion (see 
Section 4.3 i and {S/R)ies x N/b proposed by [Richter et al. 


( [2006^ for broad H I absorption systems. The motivation for 
the later being that what matters to confidently detect a broad 
absorption line is both S/N and the optical depth at the line 
centres, tq oc N/b, and so it will not be appropriate to use the 


We note that the expected numher of clusters per random sightline 
is ~ 32 (see Section j^A) , and so 162 144/32 ^ 5000 should cover 
a significant fraction of the cluster catalog in our chosen redMaPPer 
suhvolumes. 



log((5//V>,,, X N/b) 


Figure Bl. Comparison between two different methods for estimat¬ 
ing the significance level of absorption line features. The a:-axis cor¬ 
responds to the log((S/W}res X N/b) values proposed by |Richt^ 
|et al.| l [2Q06) (at least for broad H I lines) where N is in units of cm“^ 
and b in km s~^, while the y-axis corresponds to our Wr/SWr (see 
Section [43) . Blue triangles correspond to narrow H I absorption com¬ 
ponents (6 < 50 km s~^), while red pentagons con'espond to broad 
H I absorption components (6 ^ 50 km s“^), both found across the 
full Q1410 spectrum. White crosses mark components having multiple 
transitions, hence reliable (label ‘a’; see Section [43) . The grey shaded 
area con'espond to values Wr / 5Wr < 3, while the vertical dashed line 
correspond to a value of (S/N)res x N/b = 2.5 X cm-2 (km 
s“'^)“^.See AppendixjBjfor further details. 


commonly adopted formalism based on a minimum equivalent 
width threshold for unresolved lines. 

Figure |BT] shows a comparison between these two quantities 
{S/R}reE X N/b in the a;-axis and Wr/SWr in the t/-axis), for 
our sample of H I absorption lines. Blue triangles correspond to 
H I absorption components with 6 < 50 km s“^ (narrow), while 
red pentagons correspond to H I absorption components with 
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b ^ 50 km s“^ (broad) over the full Q1410 spectrum. White 
crosses mark components having multiple transitions, which we 
always account as reliable (label ‘a’; see Section [43| >. The grey 
shaded area corresponds to values Wr /<5144 < 3. 

Restricting ourselves to broad H I lines (red pentagons in 
Figure [BT| >, we see a clear correlation between these two criteria. 
The dashed vertical line in Figure [BT] corresponds to a value of 
(S/R) res X N/b ^ 2.5 X 10^^ crn”^ (km s“^)“^ for our COS 
data, which is needed to make both approaches roughly consistent 
with each other for broad H I lines. Such a value is similar to the 
one reported by Richter et al. 1 2006 ', {S/R)reB x N/b 3 x 10^^ 
cm“^ (km in their STIS data. 

For narrow Ft I lines (blue triangles in Figure |BTJ how¬ 
ever, no clear correlation is present between these two methods, 
where several components having Wr / SWr < 3, can still have 
large values of {S/R)res x N/b. We note that many of the large 
(S/R) res X N/b components are reliable (white crosses in Fig- 
ure |Bl| >, but many others at the lower end are not. In terms of 
Wr/5Wr, we see that some reliable lines fall below the limit 
of Wr/5Wr < 3. This is due to the fact that our adopted con¬ 
servative method for estimating 144 (see Section [4(3) l will give 
unrealistically large uncertainties, 5144, for unconstrained or 
saturated components. Indeed the vast majority of reliable com¬ 
ponents with 144/5144 < 3 are saturated lines for which the 
column densities are not well determined. 

Our proposed significance estimation based on 144/5144 
(and the presence of multiple transitions in a given component; 
see Section [43] ) can be straightforwardly applicable to any ab¬ 
sorption line irrespective of its Doppler parameter and ionic 
transition. Therefore it has the advantage of allowing an homoge¬ 
neous analysis. 


APPENDIX C: ESTIMATION OF REDSHIFT PATHS 
AND ABSORPTION DISTANCES 

We estimate Az{Ad, Av, 144™'“) and AX {Ad, Av, 444"““) as 
follows. First, for a given transition we just considered re¬ 
gions in the Q1410 spectrum probing rest-frame velocity dif¬ 
ferences within Av from any cluster-pair within Ad from 
the Q1410 sightline, and within our redshift range of 0.1 ^ 
« ^ 0.5. We also masked out regions with fluxes below 50% 
the value of the continuum fit, to account for the fact that 
we are usually biased against finding absorption systems on 
top of strong absorption lines. We also masked out regions 
within ±200 km s“^ from strong Galactic absorption could exist, 
namely C II A1334.53, N V AA1238.82,1242.80, O I A1302.17, 
O I* A1304.86, O I** A1306.03, Si ii AA1260.42,1304.37, 
P III A1334.81, S II AA1253.81,1259.52 and Fe ii A1260.53. 
We then calculated the minimum rest-frame equivalent width 
to observe a transition at rest-frame wavelenght, Ao, along the 
spectrum as. 


2 f'“’‘(Ad, Av, 144"“"), in the *-th spectral chunk. We then com¬ 
puted the redshift path as. 


Az{Ad, Av, wr") = X](«r“ - , 

i 

and the absorption distance as, 

max 2 

AX {Ad, Av, 44^r") = ^ ^ 


(C2) 


■\/Dni(l ± z)^ ± Da 


dz , 


(C3) 


conforming our adopted cosmology. 


APPENDIX D: EQUIVALENT WIDTHS 
DISTRIBUTIONS 

In this section we explore how do the equivalent widths dis¬ 
tributions from our sample of lines associated with cluster- 
pairs (i.e. adopting our fiducial values of Ad = 3Mpc and 
Av = 1000 km s“^), compare to those of the field expectation. 

Figure [dT| shows the redshift number density, dN/dz, of 
our different samples of absorption lines as a function of min¬ 
imum rest-frame equivalent width, 444"“", associated with our 
cluster-pairs: total H I (top left panel; brown circles), O VI (top 
right panel; green squares), narrow H I (bottom left panel; blue 
triangles) and broad H I (bottom right panel; red pentagons). The 
field distributions estimated from our Q1410 sightline are repre¬ 
sented by the dashed line with its Ala uncertainty represented 
by the darker grey region. Such field estimation comes from lines 
not being in the sample associated with cluster pairs, and its 
redshift path is estimated by excluding that of cluster-pairs too. 

We observe that the equivalent widths distributions are not 
remarkably different between lines associated with cluster-pairs 
and the field expectation, and hence the relative excess in dN/dz 
seems to remain somewhat constant as a function of 444’“'“ • Ths 
Kolmogorov-Smimov test (K-S) for the unbinned 444 distribu¬ 
tions between our cluster-pair and field samples give no statis¬ 
tically significant differences between them, but we note any 
possible real difference would be hard to detect because the 
samples are quite small (< 10). 

We conclude that, in terms of equivalent widths, our cluster- 
pair absorption line samples and the field expectations are not 
significantly different, at least from what can be obtained in our 
limited sample. 


APPENDIX E: QI410 HST/COS SPECTRUM 

Figure [ET] shows the reduced Q1410 spectrum (black line), its 
corresponding uncertainty (green lines) and our adopted pseudo¬ 
continuum fit (blue dotted line). The figure also shows our Voigt 
profile fit solutions and residuals (red lines and grey dots re¬ 
spectively; see Section [4~2l ), and their corresponding IDs and 
reliability labels (see Section[43]and Table|^. 




^ o N/R 
(S/R) 


(Cl) 


where R is the resolution of the spectrograph (taken to be 
R — 20 000), and (S/R) is the signal-to-noise ratio of the spec¬ 
trum smoothed over 2 pixels using a box-car filter (i.e. ~ 1 resolu¬ 
tion element). We then identified chunks of the spectrum satisfy¬ 
ing the criteria of 444 “a" ^ 444"““, and recorded each correspond¬ 
ing minimum and maximum redshift, 2 ™'“(Aci, Av, 444’“'") and 


APPENDIX F: VISUAL ASSOCIATION OF 
ABSORPTION LINE SYSTEMS WITH CLUSTER-PAIRS 

Figure [FT] shows an schematic view of the association between 
H I and O VI absorption line systems and our known cluster-pairs 
obtained from the redMaPPer catalog (see Section]^. These 
plots are for illustrative purposes only; in cases with multiple 
grouped cluster-pairs (see Section[±3|l the chosen redshifts for 
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Figure Dl. Redshift number density of our different samples of absorption lines associated to our cluster-pairs (i.e. using our fiducial values of 
Ad = 3 Mpc and An = 1000 km s“^), as a function of minimum rest-frame equivalent width VF™'". The top panels show the results for our total 
H I (left panel) and O VI (right panel), while the bottom panels show the results for narrow (NLA; b < 50km s“^) and broad (BLA; b ^ 50km s“^) 
H I, using the same colour/symbol convention as in Figure [^to |10| Note that bins are not independent from each other, as emphasized by the light 
coloured areas. The field distributions estimated from our Q1410 sightline are represented by the dashed line with its itlcr uncertainty represented by the 
darker grey region. The total number of lines per bin are given in grey numbers on top and bottom of each panel, for our cluster-pairs and field samples 
respectively. See Appendix|^for further details. 


defining rest-frames are arbitrary, but such choice is not used in 
our scientific analyses. 


APPENDIX G: SUMMARY TABLES FOR OUR 
INCIDENCE MEASUREMENTS 


G4 show a summary of our main dN/dz and 


Tables to 

dN/dX calculations for our different samples of absorption 
lines. 


This paper has been typeset from a TgX/ KTj^ file prepared by 
the author. 
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1245 


1250 


1255 


1260 



Observed wavelength (A) 

Figure El. Observed HST/COS FUV spectrum of QSO Q1410 (black line), and its respective uncertainty (green lines). The red line correspond to our 
fit of the spectrum (see Section [42) , while the blue dotted line corresponds to the assumed unabsorbed pseudo-continuum (i.e. including broad emission 
lines and the Milky Way’s DLA; see Section [Z^ . Vertical tick lines indicate individual absorption lines, where the numbers coiTespond to the IDs given 
in Tablej^and the letters indicate their reliability (see Section [43) . Grey points show the difference between our model and the observed data. We see 
that distribution of these residuals are consistent with the uncertainty of the data. 
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Figure El - continued 
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Figure El - continued 
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Figure El - continued 
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Figure El - continued 


© 2015 RAS, MNRAS 0Q0,[T|[24l 






























































































Towards detecting the WHIM in inter-cluster filaments 31 



s 



z=0.1588 

1000 
800 
600 
400 
200 
0 

-200 
-400 
-600 
-800 
-1000 


z=0.1760 

1000 
800 
600 
400 
200 
0 



Au (km s ^) 


s 


*0 

o 

N 



-200 

-400 

-600 

-800 

-1000 


cd 

a 

o 




HI 1216 




OVI 1032, 




OVI 1038, 








-1000 -500 


500 1000 


At; (km s ^) 



z=0.3421 

1000 
800 
600 
400 
200 
0 

-200 
-400 
-600 
-800 
I -1000 


s 



Av (km s ^) 


Figure FI. Schematic view of the clusters and cluster-pairs around the Q1410 sightline at different given redshifts (see top label in subpanels). The 
left panels correspond to the projected in the sky distribution of clusters (circles) and cluster-pairs (dashed lines) around the Q1410 sightline (yellow 
star) in co-moving Mpc. Numbers close to the circles correspond to the cluster IDs as given in Table[^ The colour of the circles correspond to the 
rest-frame velocity difference of each cluster with respect to the given redshift, according to the colour bar scale on the right of the panel. The right panels 
correspond to portions of the normalized Q1410 HST/COS spectrum within a rest-frame velocity window of ±1000 km s~^ from the given redshift, 
for the H I Lya transition (top) and the O VI AA1032,1038 Atransitions (middle and bottom, respectively). The spectrum itself is represented by the 
black lines and the uncertainty is represented by the green line. Our fits to associated absorption line systems are represented by the thick coloured lines 
according to the colour bar scale on the left of the panel, while the fits for unassociated absorption line systems (i.e. at different redshifts) are represented 
by the red lines. Absorption line ID and reliability flags are given for associated absorption line systems. 
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Towards detecting the WHIM in inter-cluster filaments 33 


Table Gl. Summary of relevant quantities for our sample of total H I absorption lines. 
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(1) Maximum transverse separation between cluster-pair axes and the Q1410 sightline (see Section lTT] . (2) Maximum velocity difference to any cluster-pair within 
Ad from the Q1410 sightline (see Section |TT] . (3) Minimum rest-frame equivalent width (see Section lsT} . (4) Redshift path (see Section [TT] and Appendix^. (5) 
Absorption distance (see Appendixj^. (6) Total number of absorption components having ^ within Ad and Av from clus ter-p airs (see Section [^^ . (7) 

Redshift number density associated to absorption components having ^ within Ad and Av from cluster-pairs (see Section 5. l|. (8) Absorption distance 

number density to absorption components having Wj: ^ within Ad and Av from cluster-pairs (analogous to dNfdz but see also Appendixj^ (9) Field 

expectation for the redshift number density (see Section [T^ . (10) Field expectation for the absorption distance number density (see Sectionj^j. (11) Excess of the 
measured redshift number density compared to the field expectation (see Sectionj^j. (12) Excess of the measured absorption distance number density compared to the 
field expectation (see Sectionj^j. 
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Table G2. Summary of relevant quantities for our sample of total O VI absorption lines. 
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q 0+4.5 

d-O —3 2 

3 7 + 4-9 
6 -^- 3 .1 

3.0 

1000 

0.07 

0.05 

0.08 

2 

aotls 

25tf, 

15+13 

lots 

9 5 + 3,8 
2-^-2,3 

9 5 + 3.6 
2-6-2.3 

3.0 

1000 

0.08 

0.05 

0.08 

2 


25 + 94 

zO-17 

15+13 

lot® 

9 5 + 3,8 
*^•^-2,3 

9 5 + 3.6 

■^■6-2.3 

3.0 

1000 

0.09 

0.05 

0.08 

2 


25 + 94 

^0-17 

15+13 

lot® 

9 5 + 3.8 
^•^-2,3 

9 5 + 3.6 

■^■6-2.3 

3.0 

1000 

0.10 

0.05 

0.08 

2 


25 + 94 

^0-17 

15+13 

lot® 

9 5 + 3.8 
^•^-2,3 

9 5 + 3.6 

■^■6-2.3 

3.0 

1000 

0.11 

0.05 

0.08 

1 


13 + 29 
13-11 

15+)3 

lot® 

1 3 + 6-1 

1 9 + 3.0 

l-‘^- 1.2 

3.0 

1000 

0.12 

0.05 

0.08 

1 

19+45 

iy_17 

13 + 29 
13_ii 

llt^^ 

0-1-8 

^-4 

1 7+^'2 

1 6 + 4-0 
1-6-1,6 

3.0 

1000 

0.14 

0.05 

0.08 

1 

19+45 

^^-17 

13 + 29 
13_ii 

8 tl 0 

5 + 7 
^-3 

2.5t®3l 

9 5 + 6.2 

•^■6-2.5 

3.0 

1000 

0.16 

0.05 

0.08 

1 

19+45 

^^-17 

13 + 29 
13_ii 

Q+10 

^-5 

5 + 7 

^-3 


2 . 5 +IJ 

3.0 

1000 

0.18 

0.05 

0.08 

1 

19+45 

iy_17 

13 + 29 
13_ii 

4 + 9 
^-3 

q-|-6 

S.lt^sV 

4 9 + 13-8 
^•^ — 4.9 


(1) Maximum transverse separation between cluster-pair axes and the Q1410 sightline (see Section lTT] . (2) Maximum velocity difference to any cluster-pair within 
Ad from the Q1410 sightline (see Section |TT] . (3) Minimum rest-frame equivalent width (see Section lsT} . (4) Redshift path (see Section [TT] and Appendix^. (5) 
Absorption distance (see Appendixj^. (6) Total number of absorption components having ^ within Ad and Av from clus ter-p airs (see Section [^^ . (7) 

Redshift number density associated to absorption components having ^ within Ad and Av from cluster-pairs (see Section 5. l|. (8) Absorption distance 

number density to absorption components having Wj- ^ within Ad and Av from cluster-pairs (analogous to dNfdz but see also Appendixj^ (9) Field 

expectation for the redshift number density (see Section [T^ . (10) Field expectation for the absorption distance number density (see Sectionj^j. (11) Excess of the 
measured redshift number density compared to the field expectation (see Sectionj^j. (12) Excess of the measured absorption distance number density compared to the 
field expectation (see Sectionj^j. 
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Table G3. Summary of relevant quantities for our sample of narrow H I (6 < 50 km s absorption lines. 


Ad 

An 

^min 

Az 

AX 

N 

dN 

dz 

dN 

dX 

dN 1 

dz Ifield 

dN 1 

dX Ifield 

dN / dN 1 
dz ' dz Ifield 

dN / dN 1 
dX ' dX Ifield 

(Mpc) 

(km s “ ^) 

(A) 










( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

3.0 

500 

0.04 

0.03 

0.05 

7 

228^1^ 

isilfe 

106 + 29 
-‘-^°-20 

75 + ^9 

'^-14 

2 2 + ^2 
2 -"^- 0,8 

2 0 +^'^ 
2 - 0 _o .8 

3.0 

630 

0.04 

0.04 

0.06 

7 


1941" 

106+29 

-‘-UD_ 2 o 

75 + ^9 

'0-14 

1 -^- 0 .7 

4 7 + 0.9 
-‘-■'- 0.6 

3.0 

794 

0.04 

0.05 

0.07 

8 


116+^^ 

^-^“-40 

106+29 

-‘-UD_ 2 o 

75 + ^9 

'0-14 

1 6+^'9 
-^■'^- 0.6 

1 6+9-9 
-‘-■“- 0.6 

3.0 

1000 

0.04 

0.06 

0.08 

10 

179+’^’^ 

120 I", 

106 + 29 
-‘-^°-20 

75 + ^9 

'0-14 

1 7+0.7 
-^■'- 0,6 

1 6+9-7 

-‘-■“-0.5 

3.0 

1313 

0.04 

0.07 

0.11 

11 

154+62 

lOSlfi 

106 + 29 
-‘-^°-20 

75 + ^9 

'0-14 

1 5 + 9'9 

1 4 + 0.6 
-‘-■^-0.4 

3.0 

1724 

0.04 

0.09 

0.13 

11 


82133 

106111 

75 + ^9 

'0-14 

T 9 + 0,5 
1 -^- 0 .4 

4 4 + 0.5 

-‘-■-‘--0.3 

3.0 

2264 

0.04 

0.12 

0.17 

11 

95 + 38 

64+26 

^^-19 

106111 

75 + ^9 

'0-14 

0 9+9-4 
'^■^-0.3 

0 9+9-4 

“■‘^-0.3 

3.0 

2972 

0.04 

0.14 

0.21 

14 

101 +^^ 

iui_27 

A« + 23 

68 _i 8 

106 + 29 
-‘-^°-20 

75 + ^9 

'0-14 

1 0+9'9 

-^■’^-0.3 

0 9+0.3 
“■‘=^-0.3 

3.0 

3903 

0.04 

0.16 

0.24 

17 

104+22 

70+21 

'^-17 

106+29 

-‘-UD_ 2 o 

75 + ^9 

'0-14 

1 0+9-9 

-^■^-0.3 

0 9+0-3 

U.y -0 2 

3.0 

5125 

0.04 

0.18 

0.27 

18 

99 + 29 

^^ — 23 

67+29 

°'-16 

106+29 

-‘-UD_ 2 o 

75 + ^9 

'0-14 

0 9+3-3 

u.y -0 2 

o.9l«;3 

3.0 

6729 

0.04 

0.19 

0.28 

18 

qq + 28 

63 + ^9 

^^-15 

106 + 29 
-‘-^°-20 

75 + ^9 

'0-14 

0 g+3’3 

u.y -0 2 

0 . 8 t«-J 

3.0 

8835 

0.04 

0.22 

0.32 

24 

110 + 2 ’^ 

76119 

106 + 29 
-‘-^°-20 

75 + ^9 

'0-14 

1 0+9-9 
-^■^- 0,2 

1 0+9-9 
-‘-■“- 0.2 

3.0 

11601 

0.04 

0.25 

0.36 

29 

115+26 

liD_ 2 i 

soils 

106111 

75 + ^9 

'0-14 

1 4 + 0,3 
-^■-^- 0.2 

4 4 + 0.3 
-‘-■-‘--0.2 

3.0 

15232 

0.04 

0.30 

0.42 

33 

112 + 2 ^ 

79 + 16 
'^_14 

loelll 

75 + ^9 

'0-14 

1 1 + 9-2 

1 1+9-2 
-^■-^- 0.2 

3.0 

20000 

0.04 

0.33 

0.47 

37 

iiaffs 

79 + 15 

'^^-13 

106 + 29 
-‘-^°-20 

75 + ^9 

'0-14 

1 1 + 9-2 
-^■-^- 0.2 

1 1+9-2 
-‘-■-‘-- 0.2 

1.0 

1000 

0.04 

0.01 

0.02 

1 

68 llf 


106 + 29 
-‘-^°-20 

75 + ^9 

'0-14 

0 6 +^'^ 
l'-°- 0.6 

0 6+^'9 

“■“-0.5 

1.4 

1000 

0.04 

0.03 

0.04 

2 

77 + 102 
‘ -51 

501®33 

106 + 29 
-‘-^°-20 

75 + ^9 

'0-14 

0 7+^'9 
l'-'-0,5 

0 7+9-9 

“■ ' -0.4 

2.1 

1000 

0.04 

0.04 

0.06 

3 

691^1 

46+45 

106+29 

-‘-UO_ 2 o 

75 + ^9 

'0-14 

0 7+9-6 
- O.A 

0 6 + 0-6 
“■“-0.4 

3.0 

1000 

0.04 

0.06 

0.08 

10 

179^^^ 

120 III 

106 + 29 

-‘-UD_ 2 o 

75 + ^9 

'0-14 

4 7 + 0,7 
- 0.6 

1 6+9-7 

-‘-■“-0.5 

4.1 

1000 

0.04 

0.06 

0.09 

10 


117l?s 

106111 

75 + ^9 

'0-14 

4 7 + 0,7 
-^■'-0.5 

1 6+9-7 

-‘-■“-0.5 

5.7 

1000 

0.04 

0.06 

0.09 

10 

174+74 

'^-54 

liolfe 

106111 

75 + ^9 

'0-14 

1 6+9-'^ 

-^■”-0.5 

1 6+9-7 

-^•“-0.5 

7.8 

1000 

0.04 

0.07 

0.10 

10 

151+66 

-L01_47 

1031" 

106+29 

-‘-UD_ 2 o 

75 + ^9 

'0-14 

1 4 + 0.6 
-^■^-0.5 

1 4 + 0-6 
-‘-■^-0.5 

10.7 

1000 

0.04 

0.08 

0.12 

11 

141 + 57 
^•^-L_42 

95111 

106 + 29 

75 + ^9 

'0-14 

T 0 + O .6 

-‘-■9-0,4 

1 0+0,5 
-‘-■“-0.4 

14.8 

1000 

0.04 

0.09 

0.14 

12 

129+^9 

izy_37 

sslll 

106 + 29 

75 + ^9 

'0-14 

T 9 + 0.5 
-1-2_0,4 

1 9+0.5 

-‘-■2-0.4 

20.3 

1000 

0.04 

0.11 

0.16 

14 

129t« 

Q O+30 
®°-23 

106 + 29 

-‘-UO_ 2 o 

75 + ^9 

'0-14 

1 2 + 9-4 
-‘-■2-0.3 

1 2+9-4 
-‘-■2-0.3 

27.9 

1000 

0.04 

0.12 

0.17 

16 

135fg 

93111 

106+29 

-‘-UD_2o 

75 + ^9 

'0-14 

1 3 + 9-4 
-‘-■^-0.3 

1 2+9-4 
-‘-■2-0.3 

38.4 

1000 

0.04 

0.15 

0.22 

19 

127+36 

-29 

87+25 

-20 

106111 

75 + ^9 

'0-14 

1 2 + 9-4 
-‘-■2-0.3 

1 2+9-4 
-‘-■2-0.3 

52.9 

1000 

0.04 

0.17 

0.24 

19 


70+22 

^®-18 

106111 

75 + ^9 

'0-14 

4 4 + 0,3 
-^■-^-0.3 

4 4 + 0.3 
-^•■^-0.3 

72.7 

1000 

0.04 

0.19 

0.28 

21 


76+21 

'°-16 

106+29 

-‘-UO_2o 

75 + ^9 

'0-14 

1 0+9-9 
-‘-■'-^- 0.2 

1 0+9-9 
-‘-■“- 0.2 

100.0 

1000 

0.04 

0.24 

0.35 

29 

119122 

84119 

106 + 29 

75lil 

4 4 + 0.3 
-‘-■-‘--0,2 

4 4 + 0.3 
-‘-■-‘--0.2 

3.0 

1000 

0.03 

0.06 

0.09 

10 

1761^ 

118111 

122 III 

87ti« 

1 5 + 9-7 
■^■‘^-0.5 

1 4 + 0-6 
-‘-■^-0.5 

3.0 

1000 

0.04 

0.06 

0.09 

10 

i^elK 

118111 

122 III 

Q> 7 "hl 8 

®'-15 

1 5 + 9-7 

■^■‘^-0.5 

1 4 + 0.6 
-‘-■^-0.5 

3.0 

1000 

0.04 

0.06 

0.09 

10 

i^elK 

118111 

11«+24 

84i;i 

1 5 + 9-7 

■^■‘^-0,5 

1 4+9-7 

-‘-■^-0.5 

3.0 

1000 

0.05 

0.06 

0.09 

10 

1761^ 

118111 

114+24 

-L- 1^-20 

82lil 

1 5 + 9-7 
-‘-■^- 0.6 

1 4+9-7 

-‘-■^-0.5 

3.0 

1000 

0.05 

0.06 

0.09 

10 

1761^, 

118111 

107+29 

-19 

77+17 

' ' -14 

1 6 + 9-9 
-^■^- 0.6 

1 5+9-7 
-‘■•‘^- 0.6 

3.0 

1000 

0.06 

0.06 

0.09 

9 

1591^2 

106111 

107111 

77+17 
' ' _14 

1 5 + 9-7 
-‘-■^- 0.6 

1 4+9-7 

-‘-■^-0.5 

3.0 

1000 

0.07 

0.06 

0.09 

9 

1591^1 

106111 

104111 

74+16 

'^-14 

1 5 + 9-3 
-‘-■^- 0.6 

1 4+9-7 
-‘-■^- 0.6 

3.0 

1000 

0.08 

0.06 

0.09 

9 

1591^^ 

106111 

104+29 

74+16 

'^-14 

1 5 + 9-3 
-‘-■^- 0.6 

1 4+9-7 
-‘-■^- 0.6 

3.0 

1000 

0.09 

0.06 

0.09 

8 

141 + 70 

94+47 

^^-33 

96lfs 

69+^9 

Dy-13 

1 5 + 9-3 
-‘-■^- 0.6 

1 4+9-7 
-‘-■^- 0.6 

3.0 

1000 

0.10 

0.06 

0.09 

7 

123+67 

143_46 

82111 

82111 

59 + 15 
09-12 

1 5 + 9-0 

-‘-■‘^-0,7 

1 4 + 0.8 
-‘-■^- 0.6 

3.0 

1000 

0.11 

0.06 

0.09 

7 

123+67 

J-4cJ_46 

82111 

79+21 

'^-17 

56 + f 5 
00-12 

1 6+9-9 

-‘-■°-0.7 

1 5 + 9-0 
-‘-■“- 0.6 

3.0 

1000 

0.12 

0.06 

0.09 

7 

123+67 

J-4cJ_46 

82111 

nlll 

51 + ^^ 
01-11 

1 7+^'9 
- 0.8 

1 6+^9 
-‘-■“-0.7 

3.0 

1000 

0.14 

0.06 

0.09 

7 

123+67 

14c5_46 

82111 

68111 

48+14 

^ 0-11 

1 R+ 1'1 
-‘-■®- 0,8 

4 7 + 1.0 
- 0.8 

3.0 

1000 

0.16 

0.06 

0.09 

7 

123+67 

14CJ_4g 

82111 

57+^9 

-14 

41 + 19 
^-^-10 

9 9 + 1.3 
2-2_i.o 

2 0+^-2 
2-“-0.9 

3.0 

1000 

0.18 

0.06 

0.09 

6 

106 I 42 

21 IS 

46+^^ 

^®-13 

33^13 

9 0 + 1.5 
z.o_i 2 

2 1 + 4-4 
^■-‘-- 1.1 


(1) Maximum transverse separation between cluster-pair axes and the Q1410 sightline (see Section [TT] . (2) Maximum velocity difference to any cluster-pair within 
Ad from the Q1410 sightline (see Section |TT] . (3) Minimum rest-frame equivalent width (see Section |5'7} . (4) Redshift path (see Section [5~T] and AppendixO. (5) 
Absorption distance (see Appendixj^. (6) Total number of absorption components having ^ within Ad and A.v from clus ter-p airs (see Section [^^ . (7) 

Redshift number density associated to absorption components having ^ within Ad and An from cluster-pairs (see Section 5. l|. (8) Absorption distance 

number density to absorption components having ^ within Ad and An from cluster-pairs (analogous to dN/dz but see also Appendixj^ (9) Field 

expectation for the redshift number density (see Sectionj^j. (10) Field expectation for the absorption distance number density (see Sectionj^j. (11) Excess of the 
measured redshift number density compared to the field expectation (see Sectionj^j. (12) Excess of the measured absoi-ption distance number density compared to the 
field expectation (see Sectionj^j. 
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Table G4. Summary of relevant quantities for our sample of broad H I (6 ^ 50 km s absorption lines. 


Ad 

An 

^min 

Az 

AX 

N 

dN 

dz 

dN 

dX 

dN 1 

dz Ifield 

dN 1 

dX Ifield 

dN / dN 1 
dz ' dz Ifield 

dN / dN 1 
dX ' dX Ifield 

(Mpc) 

(km s “ ^) 

(A) 










(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

3.0 

500 

0.04 

0.03 

0.05 

4 


86_42 

1«+13 

12+® 

7 4+0-4 
'■^-4,4 

7 0+6.0 

' •'^-4.1 

3.0 

630 

0.04 

0.04 

0.06 

4 

106+41 

7+1=4 

1«+13 

^^-10 

12+® 

6 l+^'2 
°-^-3.6 

5 7+^-^ 

' -3.4 

3.0 

794 

0.04 

0.05 

0.07 

6 


87tS 

1«+13 

^^-10 

124 

7 5 + ^' ^ 
'■^-3.9 

7 0+4-8 
' •'•'-3.7 

3.0 

1000 

0.04 

0.06 

0.08 

6 

108+^ 

7+?2 

.0+13 

^^-10 

124 

®-^-3.2 

r- 0 + 4,0 
^•^-3.0 

3.0 

1313 

0.04 

0.07 

0.11 

6 


6642^ 

.0+13 

^^-10 

12+® 

4 8+^'2 
^■® —2,5 

4 5+3'^ 
^•3-2.4 

3.0 

1724 

0.04 

0.09 

0.13 

6 

66"^ 

44+27 

^^-18 

.0+13 

^^-10 

12+? 

0 c + 2.6 

d.»-2 0 

3-6l?:=9 

3.0 

2264 

0.04 

0.12 

0.17 

7 

61 + 33 
"^-23 

41 + 22 
^-^-15 

.0+13 

^*-10 

124 

3 5 + 2-^ 

d.O-i 7 

3-3l?;= 

3.0 

2972 

0.04 

0.14 

0.21 

7 

50 + 27 

34 + 18 

.0+13 

^^-10 

12+® 


2 7+1.7 
4-'-1.4 

3.0 

3903 

0.04 

0.16 

0.24 

8 

49 + 24 

^^-17 


.0+13 

^^-10 

124 

2'+1:3 

2 6+^-® 
^•'^-1.3 

3.0 

5125 

0.04 

0.18 

0.27 

8 

44+22 

^^-15 

3+lo 

.0+13 

^^-10 

124 

2 5+1-0 
Z.O-i 2 

2 4+1-4 
^•^-1.2 

3.0 

6729 

0.04 

0.19 

0.28 

8 

41 + 20 


.0+13 

^^-10 

124 

2 4+l'‘l 

z.o_i 1 

3.0 

8835 

0.04 

0.22 

0.32 

8 

Sy-l-lS 

-13 

2+9= 

.0+13 

^^-10 

124 

2 1 + 1.3 

2 0+^'2 
^•'•’-1.0 

3.0 

11601 

0.04 

0.25 

0.36 

9 

36 + ^® 
cSD_i2 

264s^ 

.0+13 

^^-10 

12+® 

2 0+^'2 

2 0+^'^ 
^•'•'-0.9 

3.0 

15232 

0.04 

0.30 

0.42 

10 

34+^“^ 

24+10 

.0+13 

^^-10 

124 

1 9+^'^ 
-^■^-0.9 

1 9+1-° 

-^•^-0.9 

3.0 

20000 

0.04 

0.33 

0.47 

11 

00 + 13 

2i4 

.0+13 

^*^-10 

124 

1 9+l'0 

-^■^-0.9 

1 9+1-° 

-'-•^-0.9 

1.0 

1000 

0.04 

0.01 

0.02 

2 

iqC + lSO 

93 + 124 
^'^-62 

.0+13 

^^-10 

124 

- 7+10.6 
'•'-5.7 

7 c + 10.3 
'^•^-5,6 

1.4 

1000 

0.04 

0.03 

0.04 

4 


iool“ 

.0+13 

^^-10 

124 


0 .+7,0 
^•^-4.8 

2.1 

1000 

0.04 

0.04 

0.06 

5 


yy + 53 
' ' -34 

.0+13 

^^-10 

124 


6 3+^-'^ 
°’d_3.4 

3.0 

1000 

0.04 

0.06 

0.08 

6 

lOSlfj 

72tg 

.0+13 

^^-10 

124 

6 l+^'2 

°-^-3.2 

r 0 + 4.0 
^•^-3.0 

4.1 

1000 

0.04 

0.06 

0.09 

6 

losfS 

70tfs 

.0+13 

12+® 

6 0+“^'^ 
^■'^-3.1 

5.2t=3:=o 

5.7 

1000 

0.04 

0.06 

0.09 

6 

104tS 

70lfs 

.0+13 

^^-10 

12+? 

6 0+'^'^ 
^■^—3,1 

6-6l3:= 

7.8 

1000 

0.04 

0.07 

0.10 

6 

91 + 54 
^^-36 

621=5 

.0+13 

^^-10 

124 

r 9 + 3.6 

0.2-2 7 

5 0+^'“^ 
O.U-2 6 

10.7 

1000 

0.04 

0.08 

0.12 

6 


52+31 

^^-21 

.0+13 

^^-10 

124 

4 4+3.0 
^■^-2.3 

4 2 + 2-9 
^••^-2.2 

14.8 

1000 

0.04 

0.09 

0.14 

7 

7+28 

=09+28 

^^-19 

.0+13 

^^-10 

12+® 

4 3 + 2-7 

^■d_2,2 

4 2 + 2-’^ 
^••^-2.1 

20.3 

1000 

0.04 

0.11 

0.16 

7 

64+2^ 

44+24 

^^-16 

.0+13 

^^-10 

124 

q 7+2,3 
'-1.8 

3 6+2-3 

27.9 

1000 

0.04 

0.12 

0.17 

7 

59+32 
^^ — 22 

41+22 

^-^-15 

.0+13 

^^-10 

124 

3 4+2'1 
d-^-1.7 

3 3+2-1 
d-d_1.6 

38.4 

1000 

0.04 

0.15 

0.22 

8 

c 0 + 26 

ocs_i 9 

3y+18 

-13 

.0+13 

^^-10 

12+® 

3 1 + 1 ® 
d-l-1,5 

3 n+^'3 

d.U-i 4 

52.9 

1000 

0.04 

0.17 

0.24 

8 

47+23 

-16 

33tu 

.0+13 

^*-10 

12+? 

2 7 + 1.6 

'-1,3 

2 7+^-9 

^•'-1.3 

72.7 

1000 

0.04 

0.19 

0.28 

8 

41 + 21 
^^-14 

2+lo 

.0+13 

^^-10 

124 

2 4+1-^ 
^■^-1.1 

2 3+^'“^ 

2.0_i 1 

100.0 

1000 

0.04 

0.24 

0.35 

10 

41 + 18 

^^-13 

2+9= 

i8i;=o 

124 

9 q + 1 .3 

2 . 0 _i 1 

9 q+1-3 
z.o_i 1 

3.0 

1000 

0.03 

0.06 

0.09 

6 

106^42 

yi+42 

'-^-28 


334 

5 9+^'9 
^■^ — 4.1 

5.5+tl 

3.0 

1000 

0.04 

0.06 

0.09 

6 

106^42 

71+42 

'-^-28 


134 

5 9+4'® 

^■^ — 4,1 

5.5+tl 

3.0 

1000 

0.04 

0.06 

0.09 

6 

106^42 

71+42 

'-^-28 


134 

5 9+4'® 

^■^ — 4,1 

6.5+4 

3.0 

1000 

0.05 

0.06 

0.09 

6 

106^42 

71+42 

'-^-28 


134 

^■^ — 4.1 

5.5+4 

3.0 

1000 

0.05 

0.06 

0.09 

6 


71+42 
' -^-28 


134 

5 9+4’9 
^■^—4.1 

3-644 

3.0 

1000 

0.06 

0.06 

0.09 

6 

106^42 

71+42 

'-^-28 


134 

5 9+4'® 
^■^—4,1 

6.544 

3.0 

1000 

0.07 

0.06 

0.09 

6 

106^42 

71+42 

'-^-28 


134 

5 9+‘l'® 

^■^ — 4.1 

5.5+^ 

3.0 

1000 

0.08 

0.06 

0.09 

6 

106^42 

71+42 

'-^-28 


134 

^■^ — 4.1 

6-5+4 

3.0 

1000 

0.09 

0.06 

0.09 

6 

106^42 

71+42 

'-^-28 


134 

5 9+4'® 

6.5+4 

3.0 

1000 

0.10 

0.06 

0.09 

5 

0 0 + 89 

°®-38 

=iq+40 

^^-26 

14+P 

104 

6.244 

6-3+4:t 

3.0 

1000 

0.11 

0.06 

0.09 

5 

OO + 60 
®®-38 

cq+40 

^^^-26 

11 + 10 

34 

3.244 

r-7+44 

3.0 

1000 

0.12 

0.06 

0.09 

4 

yi+56 
' ^-34 

47+37 

^'-23 


34 

« « + 7.3 
®-6-6.0 

ft 1+6.8 

^•1-5.6 

3.0 

1000 

0.14 

0.06 

0.09 

4 


47+37 

^'-23 

7t=5 

5 + ^ 

^-3 

9 9-I-12'6 

9 2+^^-’^ 

3.0 

1000 

0.16 

0.06 

0.09 

2 

3641 

24+31 

■^^-16 

4+8 

^-3 

q+6 

■^-2 

9 9+20.5 

9 2 + 19,1 

3.0 

1000 

0.18 

0.06 

0.09 

2 

3641 

24+31 

•^^-16 

4+8 

^-3 

q+6 

■^-2 

9 9+20.5 
^•==' — 9,9 

9 2 + 19,1 
^•^ — 9.2 


(1) Maximum transverse separation between cluster-pair axes and the Q1410 sightline (see Section [TT] . (2) Maximum velocity difference to any cluster-pair within 
Ad from the Q1410 sightline (see Section |TT] . (3) Minimum rest-frame equivalent width (see Section |5'7} . (4) Redshift path (see Section [5~T] and AppendixO. (5) 
Absorption distance (see Appendixj^. (6) Total number of absorption components having ^ within Ad and A.v from clus ter-p airs (see Section [^^ . (7) 

Redshift number density associated to absorption components having ^ within Ad and An from cluster-pairs (see Section 5. l|. (8) Absorption distance 

number density to absorption components having ^ within Ad and An from cluster-pairs (analogous to dN/dz but see also Appendixj^ (9) Field 

expectation for the redshift number density (see Sectionj^j. (10) Field expectation for the absorption distance number density (see Sectionj^j. (11) Excess of the 
measured redshift number density compared to the field expectation (see Sectionj^j. (12) Excess of the measured absoi-ption distance number density compared to the 
field expectation (see Sectionj^j. 
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